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Abstract 

We consider here the possibihty of quantum gravity induced violation of Lorentz 
symmetry (LV). Even if suppressed by the inverse Planck mass such LV can be tested 
by current experiments and astrophysical observations. We review the effective field 
theory approach to describing LV, the issue of naturalness, and many phenomena 
characteristic of LV. We discuss some of the current observational bounds on LV, 
focusing mostly on those from high energy astrophysics in the QED sector at order 
S/Mpianck- In this context we present a number of new results which include the 
explicit computation of rates of the most relevant LV processes, the derivation of a 
new photon decay constraint, and modification of previous constraints taking proper 
account of the helicity dependence of the LV parameters implied by effective field 
theory. 
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1 Introduction 

The discovery of Lorentz symmetry was one of the great advances in the his- 
tory of physics. This symmetry has been confirmed to ever greater precision, 
and it powerfully constrains theories in a way that has proved instrumental 
in discovering new laws of physics. Moreover the mathematical structure of 
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the Lorentz group is compellingly simple. It is natural to assume under these 
circumstances that Lorentz invariance is a symmetry of nature up to arbitrary 
boosts. Nevertheless, there are several reasons to question exact Lorentz sym- 
metry. From a purely logical point of view, the most compelling reason is that 
an infinite volume of the Lorentz group is (and will always be) experimentally 
untested since, unlike the rotation group, the Lorentz group is non-compact. 
Why should we assume that exact Lorentz invariance holds when this hypoth- 
esis cannot even in principle be tested? 

While non-compactness may be a logically compelling reason to question 
Lorentz symmetry, it is by itself not very encouraging. However, there are also 
several reasons to suspect that there will be a failure of Lorentz symmetry 
at some energy or boost. One reason is the ultraviolet divergences of quan- 
tum field theory, which are a direct consequence of the assumption that the 
spectrum of field degrees of freedom is boost invariant. Another reason comes 
from quantum gravity. Profound difficulties associated with the "problem of 
time" in quantum gravity [1,2] have suggested that an underlying preferred 
time may be necessary to make sense of this physics, and general arguments 
suggest radical departures from standard spacetime symmetries at the Planck 
scale [3]. Aside from general issues of principle, specific hints of Lorentz vio- 
lation have come from tentative calculations in various approaches to quan- 
tum gravity: string theory tensor VEVs [4], cosmologically varying moduli [5], 
spacetime foam [6], semiclassical spin-network calculations in Loop QG [7,8], 
non- commutative geometry [9,10,11,12], some brane- world backgrounds [13], 
and condensed matter analogues of "emergent gravity" [14] . 

None of the above reasons amount to a convincing argument that Lorentz sym- 
metry breaking is a feature of quantum gravity. However, taken together they 
do motivate the effort to characterize possible observable consequences of LV 
and to strengthen observational bounds. Moreover, apart from any theoretical 
motivation, significant improvement of the precision with which fundamental 
symmetries are tested is always desirable. 

The study of the possibility of Lorentz violation is not new, although it has 
recently received more attention because of both the theoretical ideas just 
mentioned and improvements in observational sensitivity and reach that allow 
even Planck suppressed Lorentz violation to be detected (see e.g. [15] for an 
extensive review) . A partial list of such "windows on quantum gravity" is 

• sidereal variation of LV couplings as the lab moves with respect to a pre- 
ferred frame or directions 

• cosmological variation of couplings 

• cumulative effects: long baseline dispersion and vacuum birefringence (e.g. of 
signals from gamma ray bursts, active galactic nuclei, pulsars, galaxies) 

• new threshold reactions (e.g. photon decay, vacuum Cerenkov effect) 
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• shifted existing threshold reactions (e.g. photon annihilation from blazars, 

GZK reaction) 

• LV induced decays not characterized by a threshold (e.g. decay of a particle 
from one hehcity to the other or photon sphtting) 

• maximum velocity (e.g. synchrotron peak from supernova remnants) 

• dynamical effects of LV background fields (e.g. gravitational coupling and 
additional wave modes) 

The possibility of interesting constraints (or observations) of LV despite Planck 
suppression arises in different ways for the different types of observations. 
In the laboratory experiments looking for sidereal variations, the enormous 
number of atoms allow variations of a resonance frequency to be measured 
extremely accurately. In the case of dispersion or birefringence, the enormous 
propagation distances would allow a tiny effect to accumulate. In the new 
or shifted threshold case, the creation of a particle with mass m would be 
strongly affected by a LV term when the momentum becomes large enough 
for this term to be comparable to the mass term in the dispersion relation. 
Finally, an upper bound to electron group velocity, even if very near the speed 
of light, can severely limit the frequency of synchrotron radiation. We shall 
discuss examples of all these phenomena. 

The purpose of the present paper is twofold. First, we aim to give an introduc- 
tory overview of some of the important issues involved in the consideration of 
Lorentz violation, including history, conceptual basis and problems, observ- 
able phenomena, and the current best constraints in certain sectors. Second, 
we present a number of new results, including computations of rates of certain 
LV processes and the derivation of a new photon decay constraint. We also 
analyze the modifications of previous constraints that are required when the 
helicity dependence of the LV parameters is properly taken into account. In 
order to make the paper most useful as an introduction to LV we have placed 
much material, including most of the technical details of the new results, in 
appendices. 

The structure of this paper is as follows. The next section gives a historical 
summary of LV research, while section 3 introduces the framework for pa- 
rameterizing LV, together with the conceptual issues this raises. In section 4 
we focus on the phenomenology of dimension 5 LV in QED, and section 5 
presents the current constraints on such LV. Constraints on other sorts of LV 
are surveyed in section 6, focusing on ultra-high energy cosmic rays, and we 
close in section 7 with a discussion of future prospects. Appendix A presents 
the analysis behind the LV synchrotron constraint, appendix B derives the LV 
threshold configuration results, and appendix C includes the computation of 
rates and thresholds for some LV processes. 
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2 A brief history of some LV resecirch 

We present at this point a brief historical overview of research related to 
Lorentz violation, mentioning some influential work but without trying to be 
complete. For a more complete review see Ref. [15]. 

The idea of cosmological variation of couphng constants goes back at least to 
Milne and Dirac beginning in the 1930's [16], and continues to be of interest 
today. This would be associated with a form of LV since the spacelike surfaces 
on which the couplings are assumed constant would define a local preferred 
frame. It has recently been stimulated both by the string theoretic expectation 
that there are moduli fields, and by controversial observational evidence for 
variation of the fine structure constant a [17]. A set of related ideas goes 
under the generic name of "varying speed of fight cosmologies" (VSL), which 
includes numerous distinct formulations [18,19]. New models and observational 
constraints continue to be discussed in the literature (see e.g. [20,21]). 

Suggestions of possible LV in particle physics go back at least to the 1960's, 
when several authors wrote on that idea [22] . The possibility of LV in a met- 
ric theory of gravity was explored beginning at least as early as the 1970's 
with work of Nordtvedt and Will [23]. Such theoretical ideas were pursued in 
the '70's and '80's notably by Nielsen and several other authors on the par- 
ticle theory side [24], and by Gasperini [25] on the gravity side. A number of 
observational limits were obtained during this period [26]. 

Towards the end of the 80's Kostelecky and Samuel [27] presented evidence for 
possible spontaneous LV in string theory, and motivated by this explored LV 
effects in gravitation. The role of Lorentz invariance in the "trans-Planckian 
puzzle" of black hole redshifts and the Hawking effect was emphasized in the 
early 90's [28]. This led to study of the Hawking effect for quantum fields 
with LV dispersion relations commenced by Unruh [29] and followed up by 
others. Early in the third millennium this line of research led to work on 
the related question of the possible imprint of trans-Planckian frequencies on 
the primordial fiuctuation spectrum [30] . Meanwhile the consequences of LV 
for particle physics were being explored using LV dispersion relations e.g. by 
Gonzalez-Mestres [31]. 

Four developments in the late nineties seem to have stimulated a surge of in- 
terest in LV research. One was a systematic extension of the standard model 
of particle physics incorporating all possible LV in the renormalizable sec- 
tor, developed by Colladay and Kostelecky [32]. That provided a framework 
for computing in effective field theory the observable consequences for many 
experiments and led to much experimental work setting limits on the LV 
parameters in the Lagrangian [33]. On the observational side, the AG ASA 
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experiment reported ultra high energy (UHE) cosmic ray events beyond the 
GZK proton cutoff [34,35]. Coleman and Glashow then suggested the possi- 
bility that LV was the culprit in the possibly missing GZK cutoff [36] ^ , and 
explored many other high energy consequences of renormalizable, isotropic LV 
leading to different hmiting speeds for different particles [39]. In the fourth 
development, it was pointed out by Amelino-Camelia et al [6] that the sharp 
high energy signals of gamma ray bursts could reveal LV photon dispersion 
suppressed by one power of energy over the mass M ~ 10~^Mp, tantalizingly 
close to the Planck mass. 



Together with the improvements in observational reach mentioned earlier, 
these developments attracted the attention of a large number of researchers 
to the subject. Shortly after Ref. [6] appeared, Gambini and PuUin [7] argued 
that semiclassical loop quantum gravity suggests just such LV. ^ Following 
this work, a very strong constraint on photon birefringence was obtained by 
Gleiser and Kozameh [42] using UV light from distant galaxies. Further stim- 
ulus came from the suggestion [43] that an LV threshold shift might explain 
the apparent under-absorption on the cosmic IR background of TeV gamma 
rays from the blazar MknSOl, however it is now believed by many that this 
anomaly goes away when a corrected IR background is used [44]. 

The extension of the effective field theory (EFT) framework to incorporate LV 
dispersion relations suppressed by the ratio -E/Mpianck was performed by Myers 
and Pospelov [45]. This allowed E/Mpianck LV to be explored in a systematic 
way. The use of EFT also imposes certain relations between the LV parameters 
for different helicities, which strengthened some prior constraints while weak- 
ening others. Using EFT a very strong constraint [46] on the possibility of a 
maximum electron speed less than the speed of light was deduced from obser- 
vations of synchrotron radiation from the Crab Nebula. However, as discussed 
here, this constraint is weakened due to the helicity and particle/anti-particle 
dependence of the LV parameters 



^ Remarkably, already in 1972 Kirzhnits and Chechin [37] explored the possibility 
that an apparent missing cutoff in the UHE cosmic ray spectrum could be explained 
by something that looks very similar to the recently proposed "doubly special rela- 
tivity" [38]. 

^ Some later work supported this notion, but the issue continues to be de- 
bated [40,41]. In any case, the dynamical aspect of the theory is not under enough 
control at this time to make any definitive statements concerning LV. 
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3 Pcirametrization of Lorentz violation 

A simple approach to a phenomenological description of LV is via deformed 
dispersion relations. This approach was adopted in much work on the subject, 
and it seems to afford a relatively theory-independent framework in which to 
explore the unknown possibilities. On the other hand, not much can really 
be predicted with confidence just based on free particle dispersion relations, 
without the use of both conservation laws and interaction dynamics. Hence 
one is led to adopt a more comprehensive LV model in order to deduce mean- 
ingful constraints. In this section we discuss these ideas in turn, ending with a 
focus on the use of effective field theory, which provides a well-motivated and 
unambiguous hypothesis that can be tested. 

3.1 Deformed dispersion relations 

If rotation invariance and analyticity around p — are assumed the dispersion 
relation for a given particle type can be written as 



where E is the energy, p is hereafter the magnitude of the three-momentum, 
and 



Since this relation is not Lorentz invariant, the frame in which it applies must 
be specified. Generally this is taken to be the average cosmological rest frame, 
i.e. the rest frame of the cosmic microwave background. ^ 

Let us introduce two mass scales, M — lO^^GeV Mpianck, the putative 
scale of quantum gravity, and fi, a particle physics mass scale. To keep mass 
dimensions explicit we factor out possibly appropriate powers of these scales, 
defining the dimensionful f^'s in terms of corresponding dimensionless param- 
eters. It might seem natural that the term with n > 3 be suppressed by 
1/M"~^, and indeed this has been assumed in many works. But following this 
pattern one would expect the n = 2 term to be unsuppressed and the n = 1 
term to be even more important. Since any LV at low energies must be small, 

^ There are attempts to interpret such deformed dispersion relations as Casimir 
invariants of a new relativity group which incorporates two invariant scales, c and 
Mpi, instead of just the speed of light like in Special Relativity (SR). These theories 
are generally called Doubly (or Deformed) Special Relativity (DSR) [38]. 



E'^ = p'^ + + A{p), 




A(p) = fjip^ + f}2P^ + ?73P^ + ViP^ H 



(2) 
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such a pattern is Tintenable. Thus either there is a symmetry or some other 
mechanism protecting the lower dimension operators from large LV, or the 
suppression of the higher dimension operators is greater than 1/M"~^. This is 
an important issue to which we return in section 3.4. 

For the moment we simply follow the observational lead and insert at least 
one inverse power of M in each term, viz. 



In characterizing the strength of a constraint we refer to the rjn without the 
tilde, so we are comparing to what might be expected from Planck-suppressed 
LV. We allow the LV parameters rji to depend on the particle type, and indeed 
it turns out that they must sometimes be different but related in certain ways 
for photon polarization states, and for particle and antiparticle states, if the 
framework of effective field theory is adopted. In an even more general setting, 
Lehnert [47] studied theoretical constraints on this type of LV and deduced 
the necessity of some of these parameter relations. 

The deformed dispersion relations are introduced for individual particles only; 
those for macroscopic objects are then inferred by addition. For example, if 
N particles with momentum p and mass m are combined, the total energy, 
momentum and mass are -Etot = NE{p), ptot = -^P, and rritot = Nm, so that 
E^^^. = pI^^ + m^^^ + N'^ A{p) . Although the Lorentz violating term can be large 
in some fixed units, its ratio with the mass and momentum squared terms in 
the dispersion relation is the same as for the individual particles. Hence, there 
is no observational confiict with standard dispersion relations for macroscopic 
objects. 

This general framework allows for superluminal propagation, and spacelike 
4-momentum relative to a fixed background metric. It has been argued [48] 
that this leads to problems with causality and stability. In the setting of a LV 
theory with a single preferred frame, however, we do not share this opinion. 
We cannot see any room for such problems to arise, as long as in the preferred 
frame the physics is guaranteed to be causal and the states all have positive 
energy. 



3.2 The need for a more complete framework 



Various different theoretical approaches to LV have been taken. Some re- 
searchers restrict attention to LV described in the framework of effective field 
theory (EFT), while others allow for effects not describable in this way, such 
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as those that might be due to stochastic fluctuations of a "space-time foam" . 
Some restrict to rotationally invariant LV, while others consider also rotational 
symmetry breaking. Both true LV as well as "deformed" Lorentz symmetry 
(in the context of so-called "doubly special relativity" [38]) have been pursued. 
Another difference in approaches is whether or not different particle types are 
assumed to have the same LV parameters. 

The rest of this article will focus on just one of these approaches, namely LV 
describable by standard EFT, assuming rotational invariance, and allowing 
distinct LV parameters for different particles. This choice derives from the 
attitude that, in exploring the possible phenomenology of new physics, it is 
useful to retain enough standard physics so that clear predictions can be made, 
and to keep the possibilities narrow enough to be meaningfully constrained. 
Furthermore, of the LV phenomena mentioned in the Introduction, only dis- 
persion and birefringence are determined solely by the kinematic dispersion 
relations. Analysis of threshold reactions obviously requires in addition an as- 
sumption of energy-momentum conservation, and to impose constraints the 
reaction rates must be known. This requires knowledge of matrix elements, 
and hence the dynamics comes into play. We therefore need a complete (at 
least at low energy) theory to properly derive constraints. Only EFT currently 
satisfies this requirement. 

This approach is not universally favored (for an example of a different approach 
see [49]). Therefore we think it is important to spell out the motivation for 
the choices we have made. First, while of course it may be that EFT is not 
adequate for describing the leading quantum gravity phenomenology effects, it 
has proven very effective and flexible in the past. It produces local energy and 
momentum conservation laws, and seems to require for its applicability just 
locality and local spacetime translation invariance above some length scale. It 
describes the standard model and general relativity (which are presumably not 
fundamental theories), a myriad of condensed matter systems at appropriate 
length and energy scales, and even string theory (as perhaps most impressively 
verified in the calculations of black hole entropy and Hawking radiation rates). 
It is true that, e.g., non-commutative geometry (NCG) can lead to EFT with 
problematic IR/UV mixing, however this more likely indicates a physically 
unacceptable feature of such NCG rather than a physical limitation of EFT. 

It is worth remarking that while we choose EFT so as to work in a complete 
and well motivated framework, in fact many constraints are actually quite 
insensitive to the specific dynamics of the theory, other than that it obeys 
energy-momentum conservation. As an example, consider photon decay. In 
ordinary Lorentz invariant physics, photon decay into an electron/positron 
pair is forbidden since two timelike four-momenta (the outgoing pair) cannot 
add up to the null four momentum of the photon. With LV the photon four 
momentum can be timelike and therefore photons above a certain energy can 
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be unstable. Once the reaction is kinematically allowed, the photon lifetime is 
extremely short 1 sec) when calculated with standard QED plus modified 
dispersion, much shorter for example than the required lifetime of 10^^ seconds 
for high energy photons that reach us from the Crab nebula. Thus we could 
tolerate huge modifications to the matrix element (the dynamics) and still 
have a photon decay rate incompatible with observation. Hence, often the 
dynamics isn't particularly important when deriving constraints. 

The assumption of rotational invariance is motivated by the idea that LV 
may arise in QG from the presence of a short distance cutoff. This suggests a 
breaking of boost invariance, with a preferred rest frame ^ but not necessar- 
ily an observable breaking of rotational invariance. Note also that it is very 
difficult to construct a theory that breaks rotation invariance while preserv- 
ing boost invariance. (For example, if a spacelike four-vector is introduced to 
break rotation invariance, the four- vector also breaks boost invariance.) Since 
a constraint on pure boost violation is, barring a conspiracy, also a constraint 
on boost plus rotation violation, it is sensible to simplify with the assumption 
of rotation invariance at this stage. The preferred frame is assumed to coincide 
with the rest frame of the CMB since the Universe provides this and no other 
candidate for a cosmic preferred frame. 

Finally why do we choose to complicate matters by allowing for different LV 
parameters for different particles? First, EFT for first order Planck suppressed 
LV (see section 3.3) requires this for different polarizations or spin states, so 
it is unavoidable in that sense. Second, we see no reason a priori to expect 
these parameters to coincide. The term "equivalence principle" has been used 
to motivate the equality of the parameters. However, in the presence of LV 
dispersion relations, free particles with different masses travel on different tra- 
jectories even if they have the same LV parameters [52,53]. Moreover, different 
particles would presumably interact differently with the spacetime microstruc- 
ture since they interact differently with themselves and with each other. For 
an explicit example see [54]. Another example of this occurs in the braneworld 
model discussed in Ref. [13], and an extreme version occurs in the proposal of 
Ref. [55] in which only certain particles feel the spacetime foam effects. (Note 
however that in this proposal the LV parameters fluctuate even for a given 
kind of particle, so EFT would not be a valid description.) 



^ See however [50] for an example where (coarse grained) boost invariance is pre- 
served in a discrete model, and [51] for a study of how discreteness may be compat- 
ible with Lorentz symmetry in a quantum setting. 
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3.3 Effective field theory and LV 



In this subsection we briefly discuss some explicit formulations of LVEFT. 
First, the (minimal) standard model extension (SME) of CoUaday and Kost- 
elecky [32] consists of the standard model of particle physics plus all Lorentz 
violating renormalizable operators (i.e. of mass dimension < 4) that can be 
written without changing the field content or violating the gauge symmetry. 
For illustration, the leading order terms in the QED sector are the dimension 
three terms 

- &aVS757"V' - \H^i><T'''i^ (4) 



and the dimension four terms 

- \k''''"'Fa,F,a + ^V(Ca6 + 4b75)7" BV, (5) 



where the dimension one coefficients Hah ^-nd dimensionless k^^"^^ Cab, and 
dab are constant tensors characterizing the LV. If we assume rotational invari- 
ance then these must all be constructed from a given unit timelikc vector u"' 
and the Minkowski metric r/ab, hence ba oc Ua, Hab = 0, k"''^'^'^ oc u^"'r]'^^^^u'^ , 
Cab o(. UaUb, and dab UaUb- Such LV is thus characterized by just four num- 
bers. 

The study of Lorentz violating EFT in the higher mass dimension sector was 
initiated by Myers and Pospelov [45]. They classified all LV dimension five 
operators that can be added to the QED Lagrangian and are quadratic in the 
same fields, rotation invariant, gauge invariant, not reducible to a combination 
of lower and/or higher dimension operators using the field equations, and 
contribute terms to the dispersion relation. Just three operators arise: 

- • d){u^Fn + 7Tl7«'"^7m(Cl + C275)(« • 5)V (6) 



where F denotes the dual of F, and Ci,2 are dimensionless parameters. The 
sign of the ^ term in (6) is opposite to that in [45], and is chosen so that 
positive helicity photons have -|-^ for a dispersion coefficient (see below) . Also 
the factor 2 in the denominator is introduced to avoid a factor of 2 in the 
dispersion relation for photons. All of the terms (6) violate CPT symmetry as 
well as Lorentz invariance. Thus if CPT were preserved, these LV operators 
would be forbidden. 
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3-4 Naturalness of small LV at low energy? 

As discussed above in subsection 3.1, if LV operators of dimension n > 4 are 
suppressed as we have imagined by 1/M"~^, LV would feed down to the lower 
dimension operators and be strong at low energies [39,45,56,57], unless there 
is a symmetry or some other mechanism that protects the lower dimension 
operators from strong LV. What symmetry (other than Lorentz invariance, of 
course!) could that possibly be? 

In the Euchdean context, a discrete subgroup of the Euchdean rotation group 
suffices to protect the operators of dimension four and less from violation of 
rotation symmetry. For example [58], consider the "kinetic" term in the EFT 
for a scalar field with hypercubic symmetry, M^^ d^(j)dy(j). The only tensor M'*'' 
with hypercubic symmetry is proportional to the Kronecker delta 5^^ , so full 
rotational invariance is an "accidental" symmetry of the kinetic operator. 

If one tries to mimic this construction on a Minkowski lattice admitting a 
discrete subgroup of the Lorentz group, one faces the problem that each point 
has an infinite number of neighbors related by the Lorentz boosts. For the 
action to share the discrete symmetry each point would have to appear in in- 
finitely many terms of the discrete action, presumably rendering the equations 
of motion meaningless. 

Another symmetry that could do the trick is three dimensional rotational 
symmetry together with a symmetry between different particle types. For ex- 
ample, rotational symmetry would imply that the kinetic term for a scalar 
field takes the form {dttpY — (^{di'PY^ for some constant c. Then, for multiple 
scalar fields, a symmetry relating the fields would imply that the constant c 
is the same for all, hence the kinetic term would be Lorentz invariant with c 
playing the role of the speed of light. Unfortunately this mechanism does not 
work in nature, since there is no symmetry relating all the physical fields. 

Perhaps under some conditions a partial symmetry could be adequate, e.g. 
grand unified gauge and/or super symmetry. In fact, some recent work [59,60] 
presents evidence that supersymmetry can indeed play this role. Supersym- 
metry (SUSY) here refers to the symmetry algebra that is a 'square root' 
of the spacetime translation group. The nature of this square root depends 
upon the Minkowski metric, so is tied to the Lorentz group, but it does not 
require Lorentz symmetry. Nibbelink and Pospelov showed in Ref. [59], us- 
ing the superfield formalism, that the SUSY and gauge symmetry preserving 
LV operators that can be added to the SUSY Standard Model first appear 
at dimension five. This solves the naturalness problem in the sense of having 
Planck-suppressed dimension five operators without accompanying huge lower 
dimension LV operators. However, it should be noted that in this scenario the 
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Oijr') terms in the particle dispersion relations are suppressed by an additional 
factor of TV? jM"^ compared to (2,3). 

In a different analysis, Jain and Ralston [60] showed that if the Wess-Zumino 
model is cut off at the energy scale A in a LV manner that preserves SUSY, 
then LV radiative corrections to the scalar particle self-energy are suppressed 
by the ratio m^/A^, unlike in the non-SUSY case [57] where they diverge 
logarithmically with A. This is another example of SUSY making approximate 
low energy Lorentz symmetry natural in the presence of large high energy LV. 

Of course SUSY is broken in the real world. LV SUSY QED with softly bro- 
ken SUSY was recently studied in [61]. It was found there that, upon SUSY 
breaking, the dimension five SUSY operators generate dimension three oper- 
ators large enough that the dimension five operators must be suppressed by a 
mass scale much greater than Mpianck- In this sense, the naturalness problem 
is not solved in this setting (although it is not as severe as with no SUSY). If 
CPT symmetry is imposed however, then all the dimension five operators are 
excluded. After soft SUSY breaking, the dimension six LV operators generate 
dimension four LV operators that are currently not experimentally excluded. 
Hence perhaps this latter scenario can solve the naturalness problem. But 
again, the SUSY LV operators considered in [61] do not give rise to the type 
of dispersion corrections we consider, hence the astrophysical bounds discussed 
here are not relevant to that SUSY model. 

At this stage we assume the existence of some realization of the Lorentz sym- 
metry breaking scheme upon which constraints are being imposed. If none 
exists, then our parametrization (3) is misleading, since there should be more 
powers of 1/M suppressing the higher dimension terms. In that case, current 
observational limits on those terms do not significantly constrain the funda- 
mental theory. 



4 Phenomenology of QED with dimension 5 Lorentz violation 

We now discuss in general terms the new phenomena arising when the extra 
dimension five operators in (6) arc added to the QED lagrangian. This lays 
the groundwork for the specific constraints discussed in the next section. 

Free particles and tree level interactions can be analyzed without specifying 

the underlying mechanism that (we assume) protects approximate low energy 
Lorentz symmetry, hence we restrict attention to these. In principle radia- 
tive loop corrections cannot be avoided, but their treatment requires some 
commitment as to the specific mechanism. 
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The appearance of higher time derivatives in the Lagrangian brings with it, 
if treated exactly, an increase in the number of degrees of freedom. In the 
EFT framework however it is natural to truncate the theory via perturbative 
reduction to the degrees of freedom that exist without the higher derivative 
terms (see e.g. [62]). Although not discussed exphcitly, it is this truncation we 
have in mind in what follows. 



4.1 Free particle states 



4.1.1 Photons 

In Lorentz gauge, d^A/j, — 0, the free field equation of motion for in the 
preferred frame from (6) is 

= ^e^p,su^u ■ dfF^\ (7) 



For a particle travelling in the z-direction, the ^0,3 equations are the usual 
□Ao,3 = and hence the residual gauge freedom can still be applied to set 
Aq — As — 0, leaving two transverse polarizations with dispersion [45] 

^l = k'± hk\ (8) 



The photon subscripts ± denote right and left polarization, e± = ± ze^, 
which have opposite dispersion corrections as a result of the CPT violation in 
the Lagrangian. 



4.1.2 Fermions 

We now solve the modified free field Dirac equation for the electron and 
positron eigenspinors and find the corresponding dispersion relations. We shall 
see that there exists a basis of energy eigenstates that also have definite he- 
licity p • J, hence helicity remains a good quantum number in the presence of 
the LV dimension five operators. 

Beginning with the electron, we assume the eigenspinor has the form Usip)e~^P'^ 
where p is the 4-momentum vector, assumed to have positive energy, and 
s = ±1 denotes positive and negative helicity. The field equation from (6) in 
the chiral basis implies 

ArUr = muL, AlUl = muR (9) 
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where 

Ar^l = £; t p • a - v±EV2M, (10) 



with 'r]± = 2(Ci ± (2) (the upper signs refer to ^4^) and a are the usual Pauh 
spin matrices. If the spinors are hehcity eigcnspinors wc can replace p • a with 
simply sp. The dispersion relation is then ArAl — m^, or 

E' = m^+P' + ^E\E + sp) + ^E\E - sp) - (11) 



Moreover, the spinors ur and ul are proportional to each other, with ratio 
ul/ur — Ar/Al. The eigenspinors for the electron can thus be taken as 



\ 



_(^E-sp---^ 



Xs{p) 



y \/AEXs{p) ) \ E + sp - ^Xs{p) 



[12) 



where Xs is the two-component eigenspinor of p ■ a. with eigenvalue s. Note 
that the dispersion relation ArA^ = implies that either Ar and Ai are 
both positive, or they are both negative. The definition (10) shows that when 
the energy E is positive and much smaller than M, at least one of the two is 
always positive, hence they are both positive. Then from (9) we see that ur 
and Ul are related by a positive real multiple, implying that the square roots 
in (12) have a common sign. 

The normahzation of this spinor is ul{p)Us{p) = Ar + Al = 2E — 2(iE'^/M. 
If the usual factor (2£')~^/^ is included in the momentum integral in the field 
operator, the spinors should be normalized to 2E to ensure the canonical com- 
mutation relations hold with the standard annihilation and creation operators 
assumed. Since we work only at energies E -C M, the correction is small and 
may be neglected. 

The constraints we shall discuss arise from processes in which the energy 
satisfies m <^ E <^ M. Then to lowest order in m and £'^/M the dispersion 
relation (11) for positive and negative helicity electrons takes the form [45] 

El=p' + w? + r^±^. (13) 



where we have replaced by p in the last term, which is valid to lowest order. 
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At E ^ m the helicity states take the approximate form 



X (p) 



(14) 



These are almost chiral, with mixing still controlled by the mass as in the 
usual Lorentz invariant case. 

To find the dispersion relation and spinor wavefunctions for positrons we can 
use the hole interpretation. A positron of energy, momentum, and spin an- 
gular momentum {E, p, J) corresponds to the absence of an electron with 
{—E, — p, — J). This electron state has the same helicity as the positron state, 
since (— p) ■ (—J) = p- J. Hence the positron dispersion relation is obtained by 
the replacement {E,s) {—E,s) in (11). This is equivalent to the replace- 
ment f]± — > —Tjzf [63], from which we conclude that the LV parameters for 
positrons are related to those for electrons by 

^positron _ _^electron ^j^g-^ 



The spinor wavefunction of the {—E, — p, —J) electron state — which is also 
the wavefunction multiplying the positron creation operator in the expansion 
of the fermion field operator — has the opposite spin state compared to the 
(£^, p,J) electron, so in place of Xs one has X-s- Since the energy —E is 
negative, both Ar and Al are now negative. Hence the positron eigenspinors 
take the form 



Vs{p) = 



' ^\X-s{p) 

,^-^\A^\X-s{p) ^ 



V - ^ E - sp + ^^x-sip) J 



(16) 



At high energy these take the approximate form 



' x-{p) ^ 

K-^X-{p)j 



V {p) 




(17) 



This completes our analysis of the free particle states. 
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4-2 LV signatures: free particles 

4.2.1 Vacuum dispersion and birefringence 

The photon dispersion relation (8) entails two free particle phenomena that 
can be used to look for Lorentz violation and to constrain the parameter ^: 

(i) The propagation speed depends upon both frequency and polarization, 
which would produce dispersion in the time of flight of different frequency 
components of a signal originating at the same event. 

(ii) Linear polarization direction is rotated through a frequency-dependent 
angle, due to different phase velocities for opposite hclicities. This vacuum 
birefringence would rotate the polarization direction of monochromatic radia- 
tion, or could depolarize linearly polarized radiation composed of a spread of 
frequencies. 

4-2.2 Limiting speed of charged particles 

The electron or positron dispersion relations (13) imply that if the r^-parameter 
is negative for a given helicity, that helicity has a maximum propagation 
speed which is strictly less than the low frequency speed of light. This LV 
phenomenon would limit the frequency of synchrotron radiation produced by 
that helicity. This is not strictly a free particle effect since it involves accel- 
eration in an external field and radiation, but the essential LV feature is the 
hmiting speed. In the Appendix section A we review the EFT analysis of this 
phenomenon. 

4.3 LV signatures: particle interactions 

There are a number of LV effects involving particle interactions that do not 
occur in ordinary Lorentz invariant QED. These effects include photon de- 
cay (7 — > e'^e~), vacuum Cerenkov radiation and helicity decay (e~ — > e~7), 
fermion pair emission (e~ e~e~e~^), and photon splitting (7 w-f). Ad- 
ditionally, the threshold for the "photon absorption" (77 e'^e~) is shifted 
away from its Lorentz invariant value. 

All but photon splitting have nonvanishing tree level amplitudes, hence can be 

considered without getting involved in the question of the UV completion of 
the theory. Once loop amplitudes are considered, the UV completion becomes 
important, and as discussed in section 3.4 can produce large effects at low 
energy unless there is fine tuning or a symmetry protection mechanism. 
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The threshold effects can occur at energies many orders of magnitude below 
the Planck scale. To see why, note that thresholds are determined by particle 
masses, hence if the 0{E/M) term in Eqs. (8) or (13) is comparable to the 
electron mass term, m^, one can expect a significant threshold shift. For LV 
coefficients of order one this occurs at the momentum 



which gives a rough idea of the energies one needs to reach in order to put 
constraints of order one on the Lorentz violations considered here. ^ 

To use the anomalous decay processes for constraining LV one needs to know 
the threshold energy (if any) and how the rate depends on the incoming parti- 
cle energy and the LV parameters. The details of computing these thresholds 
and rates are relegated to the Appendix sections B and C. Here we just cite 
the main results. 

4-3.1 Photon decay 

A photon with energy above a certain threshold can decay to an electron and a 
positron, 7 — > e'^e~ . The threshold and rate for this process generally depend 
on all three LV parameters ^ and ri± . Earlier use of the photon decay constraint 
with the dispersion relation (8) did not allow for helicity and particle/anti- 
particle dependence of the LV parameters. To obtain constraints on just two 
parameters, but consistent with EFT, we can focus on processes in which 
only either 77+ or rj- is involved, namely reactions in which the positron has 
opposite helicity to the electron. For example, if the electron has positive 
helicity then its LV parameter is 77+, and according to (15) the LV parameter 
for the negative helicity positron is —r]+. 

At threshold, the final particle momenta are all ahgned (cf. Appendix B). 
Since the incoming photon has nonzero spin, angular momentum cannot be 
conserved if the electron and positron have opposite helicities. However, the 
momenta above threshold need not be aligned, so that angular momentum 
can be conserved with opposite helicities at any energy above threshold. In 
this case, the rate vanishes at threshold and is suppressed very close to the 
threshold, but above the threshold the rate quickly begins to scale as E^/M, 
where E is the photon energy. ^ For a 10 TeV photon, this corresponds to 
a decay time of order 10~® seconds (where we have taken into account the 
extra suppression coming from overall numerical factors, see section C.2 in 

^ See also [53] for a generalization to arbitrary order of Lorentz violation and dif- 
ferent particle sectors. 

^ This corrects a prior assertion [64,53] that the rate scales as E. 




(18) 
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the Appendix). The rapidity of the reaction (in comparison to the required 
hfetime of 10^^ seconds for an observed Crab photon) imphes that a threshold 
constraint, i.e. a bound on LV coefficients such that the reaction does not 
happen at all, will be extremely close to the constraint derived by requiring 
the hfetime to be longer than the observed travel time. 

4-3.2 Vacuum Cerenkov radiation and helicity decay 

The process e~ — > e~7 can either preserve or flip the helicity of the electron. If 
the electron helicity is unchanged we call it vacuum Cerenkov radiation, while 
if the helicity changes we call it helicity decay. Vacuum Cerenkov radiation 
occurs only above a certain energy threshold, and above that energy the rate 
of energy loss scales as p^/M, which implies that a 10 TeV electron would emit 
a significant fraction of its energy in 10~^ seconds. Above the threshold energy, 
constraints derived simply from threshold analysis alone are again reliable. 

In contrast, if the positive and negative helicity LV parameters for electrons 
are unequal, say > then decay from negative to positive helicity can 
happen at any energy, i.e. there is no threshold. ^ However it can be shown (see 
Appendix C.7) that assuming generic order unity values of the LV parameters, 
the rate is extremely small until the energy is comparable to the threshold for 
Cerenkov radiation. Above this energy, the rate is ~ e^m? /p independent 
of the LV parameters, which (taking into account all the numerical factors) 
imphes that a 10 TeV electron would fiip hehcity in about 10~^ seconds. Thus, 
above this "effective threshold" , constraints on helicity decay derived from the 
value of the effective threshold alone are reliable. 



4-3.3 Fermion pair emission 

The process e~ — > e~e~e"'" is similar to vacuum Cerenkov radiation or helicity 
decay, with the final photon replaced by an electron-positron pair. This reac- 
tion has been studied previously in the context of the SME [48,65]. Various 
combinations of helicities for the different fermions can be considered indi- 
vidually. If we choose the particularly simple case (and the only one we shall 
consider here) where all electrons have the same helicity and the positron has 
the opposite helicity, then according to Eq. 15 the threshold energy will de- 
pend on only one LV parameter. In Appendix C.6 we derive the threshold for 
this reaction, finding that it is a factor ~ 2.5 times higher than that for soft 



The lack of a threshold can easily be seen by noting that the dispersion curves 
E{p) of the two helicity electrons can be connected by a null vector for any energy. 
Therefore there is always a way to conserve energy-momentum with photon emission 
(the photon four- vector is almost null even with LV). 
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vacuum Ccrcnkov radiation. The rate for the reaction is high as well, hence 
constraints may be imposed using just the value of the threshold. 

4.3.4 Photon splitting 

The photon splitting processes 7 — > 27 and 7 87, etc. do not occur in 
standard QED. Although there are corresponding Feynman diagrams (the 
triangle and box diagrams), their amplitudes vanish. In the presence of Lorentz 
violation these processes are generally allowed when ,^ > 0. However, the 
effectiveness of this reaction in providing constraints depends heavily on the 
decay rate. 

Aspects of vacuum photon splitting have been examined in [53,66]. An esti- 
mate of the rate, independent of the particular form of the Lorentz violating 
theory and neglecting the polarization dependence of the photon dispersion 
relation, was given in Ref. [53]. It was argued that a lower bound on the life- 
time is S~^E~^, where 5 is a Lorentz violating factor which for LV from the 
Myers- Pospelov lagrangian is 5 ~ ^E'/M. It was imphcitly assumed in this ar- 
gument that no large dimensionless ratios occur in the rate. However, a recent 
paper [67] shows that such a large ratio can indeed occur. 

Ref. [67] analyzes the case where there is no Lorentz violation in the electron- 
positron sector. Then the (Lorcntz-invariant) one-loop Euler-Heisenberg La- 
grangian characterizes the photon interaction, even in the "rest frame" of 
the incoming photon. Using this interaction it is argued that the lifetime is 
much shorter, by a factor {me/E^)^6~^, than the lower bound mentioned in 
the previous paragraph. The extra factor of 5 in the rate is compatible with 
the analysis of [53] since there only the minimum number of factors of S was 
determined. However, the possibile appearance of a large dimensionless factor 
like {E^/nie)^ in the rate was overlooked in [53]. 

Using 50 TeV photons from the Crab nebula, about 10^^ seconds away, the 
analysis of Ref. [53] would have implied that the constraint on ^ can be no 
stronger than ^ < 10'\ and even this is not competitive with the other con- 
straints. However, for a 50 TeV photon the extra factor in the lifetime from the 
analysis of Ref. [67] is of order 10"^'^, which would produce a bound ^ < 10^^. 
Note however that this analysis so far neglects the hehcity dependence of the 
photon dispersion, and assumes no Lorentz violation in the electron-positron 
sector. 

4-3.5 Photon absorption 

A process related to photon decay is photon absorption, 77 — > e~^e~. Unlike 
photon decay, this is allowed in Lorentz invariant QED. If one of the photons 
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has energy ujq, the threshold for the reaction occurs in a head-on coUision with 
the second photon having the momentum (cquivalently energy) ki^i = /uoq. 
For kui = 10 TeV (which is relevant for the observational constraints) the soft 
photon threshold ujq is approximately 25 meV, corresponding to a wavelength 
of 50 microns. 

In the presence of Lorentz violating dispersion relations the threshold for this 
process is in general altered, and the process can even be forbidden. Moreover, 
as noticed by Kluzniak [68] , in some cases there is an upper threshold beyond 
which the process does not occur. ^ The lower and upper thresholds for photon 
annihilation as a function of the two parameters ^ and r/, were obtained in [53], 
before the helicity dependence required by EFT was appreciated. As the soft 
photon energy is low enough that its LV can be ignored, this corresponds to 
the case where electrons and positrons have the same LV terms. The analysis 
is rather complicated. In particular it is necessary to sort out whether the 
thresholds are lower or upper ones, and whether they occur with the same or 
different pair momenta. 

The photon absorption constraint, neglecting helicity dependent effects, came 
from the fact that LV can shift the standard QED threshold for annihilation 
of multi-TeV 7-rays from nearby blazars, such as Mkn 501, with the ambient 
infrared extragalactic photons [68,70,71,53,72,73,74]. LV depresses the rate of 
absorption of one photon helicity, and increases it for the other. Although 
the polarization of the 7-rays is not measured, the possibility that one of the 
polarizations is essentially unabsorbed appears to be ruled out by the obser- 
vations which show the predicted attenuation [74] . The threshold analysis has 
not been redone to allow for the helicity and particlc/anti-particic dependence 
of 77. The motivation for doing so is not great since the constraint would at best 
not be competitive with other constraints, and the power of the constraint is 
limited by our ignorance of the source spectrum of emitted gamma rays. 



5 Constraints on LV in QED at 0{E/M) 



In this section we discuss the current observational constraints on the LV 
parameters ^ and ri± in the Myers- Pospelov extension of QED, focusing on 
those coming from high energy processes. The highest observed energies occur 
in astrophysical settings, hence it is from such observations that the strongest 
constraints derive. 

The currently most useful constraints are summarized in Fig. 1. The allowed 

^ A detailed investigation of upper thresholds was carried out by the present authors 
in [53,69]. Our results agree with those of [68] only in certain limiting cases. 
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Fig. 1. Constraints on LV in QED at 0{E/M) on a log- log plot. For negative 
parameters minus the logarithm of the absolute value is plotted, and region of 
width 10"^*^ is excised around each axis. The constraints in solid lines apply to ^ 
and both •q±, and are symmetric about both the ^ and the rj axis. At least one of 
the two pairs {ri± , ^) must lie within the union of the dashed bell-shaped region and 
its reflection about the ^ axis. The IC and synchrotron Cerenkov lines are truncated 
where they cross. 

parameter space corresponds to the dark region of the Figure. These con- 
straints strongly bound Lorentz violation at order 0{E/M) in QED, assum- 
ing the framework of effective field theory holds. While the natural magnitude 
of the photon and electron coefficients ^, r) would be of order unity if there 
is one power of suppression by the inverse Planck mass, the coefficients are 
now restricted to the region |^| < 10~^ by birefringence and \ri±\ < 10~^ by 
photon decay. The narrower bell-shaped region bounded by the dashed lines is 
determined by the combination of synchrotron and Cerenkov constraints and 
applies only to one of the four pairs {±.ri±,^). Equivalently, the union of this 
bell-shaped region with its reflection about the ^ axis applies to one of the 
two pairs {r]±,^). We shall now briefly discuss how each constraint is obtained, 
leaving the details for the Appendices. 



5. 1 Photon time of flight 



Photon time of flight constraints [75] limit differences in the arrival time at 
Earth for photons originating in a distant event [76,6]. Time of flight can vary 
with energy since the LV term in the group velocity is ^k/M. The arrival time 
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difference for wave-vectors ki and k2 is thus 



At = e(^2 - ki)d/M, 



(19) 



which is proportional to the energy difference and the distance travelled. Con- 
straints were obtained using the high energy radiation emitted by some gamma 
ray bursts (GRB) and active galaxies of the blazars class. The strength of such 
constraints is typically ^ < O(IOO) or weaker [75]. 

A possible problem with the above bounds is that in an emission event it is 
not known if the photons of different energies are produced simultaneously. 
If different energies are emitted at different times that might mask a LV sig- 
nal. One way around this is to look for correlations between time delay and 
redshift, which has been done for a set of GRB's in [77]. Since time of flight 
delay is a propagation effect that increases over time, a survey of GRB's at 
different redshifts can separate this from intrinsic source effects. This enables 
constraints to be imposed (or LV to be observed) despite uncertainty regarding 
source effects. While bounds derived in this way are, at present, weaker than 
the bound presented here, they are also more robust. One might also think 
that the source uncertainties in GRB's could be mitigated by looking at high 
energy narrow bursts, which by definition have nearly simultaneous emission. 
However, for high energy narrow bursts the number of photons per unit time 
can be very low, thereby limiting the shortest detectable time lag [78]. 

The simultaneity problem may also be avoided if one uses the EFT dispersion 
relation (8) for photons of the same energy. Then one can consider the velocity 
difference of the two polarizations at a single energy [63] 



In principle this leads to a constraint at least twice as large as the one arising 
from energy differences which is also independent of any intrinsic time lag 
between different energy photons. [It is possible, of course, that there is an 
intrinsic helicity source effect, however this seems unlikely] In Fig. 1 we use 
the EFT improvement of the constraint of Biller et al. [75], obtained using the 
blazar Markarian 421, an object whose distance is reliably known, which yields 
1^1 < 63 [63]. Note however that this bound assumes that both polarizations 
are observed. If for some reason (such as photon decay of one polarization) only 
one polarization is observed, then the bound shown in is Fig. 1 weakened by 
a factor of two. In any case the time of flight constraint remains many orders 
of magnitude weaker than the birefringence constraint so is rather irrelevant 
from an EFT perspective. 



At^2\C\k/M. 



(20) 
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5.2 Birefringence 



The birefringence constraint arises from the fact that the LV parameters for 
left and right circular polarized photons are opposite (8). The phase velocity 
thus depends on both the wavevector and the hclicity. Linear polarization is 
therefore rotated through an energy dependent angle as a signal propagates, 
which depolarizes an initially linearly polarized signal comprised of a range 
of wavevectors. Hence the observation of linearly polarized radiation coming 
from far away can constrain the magnitude of the LV parameter. 

In more detail, with the dispersion relation (8) the direction of linear polar- 
ization is rotated through the angle 



for a plane wave with wave-vector k over a propagation time t. The difference 
in rotation angles for wave-vectors ki and k2 is thus 



where we have replaced the time t by the distance d from the source to the 
detector (divided by the speed of light). Note that the effect is quadratic in 
the photon energy, and proportional to the distance travelled. The constraint 
arises from the fact that once the angle of polarization rotation differs by 
more than 7r/2 over the range of energies in a signal, the net polarization is 
suppressed. 

This effect has been used to constrain LV in the dimension three (Chern- 
Simons) [79], four [80] and five terms. The strongest currently reliable con- 
straint on the dimension five term was deduced by Gleiser and Kozameh [42] 
using UV light from distant galaxies, and is given by |^| 2 x 10~^. The much 
stronger constraint |C| ^ 2 x 10~^^ was derived [63,81] from the report [82] of 
a high degree of polarization of MeV photons from GRB021206. However, the 
data has been reanalyzed in two different studies and no statistically signifi- 
cant polarization was found [83] . 

5.3 The Crab nebula 

Apart from the two constraints just discussed, all the others reported in Fig. 1 
arise from observations of the Crab nebula. This object is the remnant of a 
supernova that was observed in 1054 A.D., and lies only about 1.9 Kpc from 
Earth. It is characterized by the most energetic QED processes observed (e.g. 



9{t) = [uj+{k) - uj-{k)] t/2 = ^kH/2M 



(21) 



(22) 
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it is the source of the highest energy gamma rays) and is very well studied. In 
contrast to the distant sources desirable for constraints based on the cumu- 
lative effects of Lorentz violation, the Crab nebula is nearby. This proximity 
facilitates the detection of the low fluxes characteristic of the emission at the 
highest energies, which is particularly useful for the remaining constraints. Be- 
fore undertaking a discussion of the constraints we first summarize the nature 
of the Crab emission. 

The Crab nebula is a bright source of radio, optical. X-ray and gamma-ray 
emission. It exhibits a broad spectrum characterized by two marked humps. 
This spectrum is consistently explained by a combination of synchrotron emis- 
sion by a high energy wind of electrons and positrons, and inverse Compton 
scattering of the synchrotron photons (plus perhaps 10% other ambient pho- 
tons) by the same charges [84,85]. No other model for the emission is under 
consideration, other than for producing some fraction of the highest energy 
photons. For the constraints discussed here we assume that this standard syn- 
chrotron/inverse Compton model is correct. In contrast to our previous work, 
here we take fully into account the role of the positrons, as well as the helicity 
and particle/anti-particle dependence of the LV parameters. This complicates 
matters and changes the nature of the constraints, weakening some aspects 
and strengthening others. 

The inverse Compton gamma ray spectrum of the Crab extends up to energies 
of at least 50 TeV. The synchrotron emission has been observed to extend at 
least up to energies of about 100 MeV [84], just before the inverse Compton 
hump begins to contribute to the spectrum. In standard Lorentz invariant 
QED, 100 MeV synchrotron radiation in a magnetic field of 0.6 mC would be 
produced by electrons (or positrons) of energy 1500 TeV. The magnetic field 
in the emission region has been estimated by several methods which agree 
on a value between 0.15-0.6 mG (see e.g. [86] and references therein). Two 
of these methods, radio synchrotron emission and equipartition of energy, are 
insensitive to Planck suppressed Lorentz violation, hence we are justified in 
adopting a value of this order for the purpose of constraining Lorentz violation. 
We use the largest value 0.6 mG for B, since it yields the weakest constraint 
for the synchrotron radiation. 



5.4 Photon decay 

The observation of 50 TeV gamma rays emitted from the Crab nebula implies 
that the threshold for photon decay for at least one helicity must be above 50 
TeV. By considering decays in which the electron and positron have opposite 
helicity we can separately constrain both (r;+,C) and (r^-,0 (cf- section 4.3.1 
and Appendix C) . The allowed region is the union of those for which positive 
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and negative helicity photons do not decay. This yields the constraint shown in 
Figure 1. The complete form of the aUowed region was determined mimericahy. 
However in the region |^| < 10""^ defined by the birefringence constraint, ^ can 
be neglected compared to r]±, in which case the constraint takes the analytic 
form |77±| < eVSrri^M/kff^. With kth = 50Tev this evaluates to |77±| < 0.2. 

5.5 Vacuum Cerenkov — Inverse Compton electrons 

The inverse Compton (IC) Cerenkov constraint uses the electrons and positrons 
of energy up to 50 TeV inferred via the observation of 50 TeV gamma rays 
from the Crab nebula which are explained by IC scattering. Since the vacuum 
Cerenkov rate is orders of magnitude higher than the IC scattering rate, that 
process must not occur for these charges [39,53] for at least one of the four 
charge species (plus or minus helicity electron or positron), and for either pho- 
ton helicity. The excluded region in parameter space is thus symmetric about 
the 77-axis, but applies only to one of the four parameters ±ri±. That is, there 
is a constraint that must be satisfied by either the pair (I77+I, |^|) or the pair 

The absence of the soft Cerenkov threshold up to 50 TeV produces the dashed 
vertical IC Cerenkov line in Fig. 1 (see Appendix C.5 for a more detailed dis- 
cussion of the threshold). One can see from (C.21) that this yields a constraint 
on f] of order (11 TeV/50 TeV)"^ ~ 10^^. By itself this imposes no constraint at 
all on the parameters r]±, since one of the two parameters ±?7_|_ always satisfies 
it, as does one of ±?7_. However, this is not the whole story. 

For parameters satisfying both ^ < —3r] and C < V ^he Cerenkov threshold 
is "hard", since it involves emission of a high energy photon [53,72]. This 
constraint depends upon both ^ and rj, and is a curve on the parameter space. 
We do not indicate this constraint in Fig. 1 since it is superseded by the 
synchrotron-hard Cerenkov constraint discussed below. It imposes a lower 
bound on one of |77±| once |^| is large enough. 

5. 6 Fermion pair emission 

If we knew that all electrons and positrons were stable at 50 TeV, then the 
threshold for all of them to undergo fermion pair emission would necessarily 
be over 50 TeV. Using the processes e" — > e~e~e'^ and e"*" — > e'^e'^e~ with 
the helicities chosen so that only one of r]± is involved in the reaction, elec- 
tron/positron stability would lead to the constraint Ir^-tl < 0.16. This is a 
factor of 16 higher than that from the soft Cerenkov threshold, and roughly 
the same as the photon decay constraint. However, at least without further 
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detailed analysis of the Crab spectrum, all we can say is that at least one 
of the fermion types is able to produce the 50 TeV IC radiation. By itself 
this imposes no constraint at all, since for each helicity either the electron or 
the positron of the opposite helicity is stable to this particular pair emission 
process. 

5.7 Synchrotron radiation 

An electron or positron with a negative value of r/ has a maximal group velocity 
less than the low energy speed of light, hence there is a maximal synchrotron 
frequency o;™^'' (A. 3) that it can produce, regardless of its energy [46] (for 
details see Appendix A). Thus for at least one electron or positron hehcity o;™^ 
must be greater than the maximum observed synchrotron emission frequency 
a;obs- This yields the constraint 



The strongest constraint is obtained in the case of a system that has the 
smallest ratio. This occurs in the Crab nebula, which emits synchrotron 

radiation up to 100 MeV and has a magnetic field no larger than 0.6 mC in the 
emitting region. Thus we infer that at least one of ±77^ must be greater than 
—7 X 10~^. This constraint is shown as a dashed line in Fig. 1. As with the soft 
Cerenkov constraint, by itself this imposes no constraint on the parameters, 
but as discussed in the next subsection it plays a role in a combined constraint 
using also the Cerenkov efi^ect. 

5.8 Vacuum Cerenkov effect — synchrotron electrons 

The existence of the synchrotron producing charges can be exploited to extend 
the vacuum Cerenkov constraint. For a given rj satisfying the synchrotron 
bound, some definite electron energy £'synch(^) must be present to produce the 
observed synchrotron radiation. (As explained in Appendix A, this is higher 
for negative rj and lower for positive rj than the Lorcntz invariant value [46].) 
Values of |(^| for which the vacuum Cerenkov threshold is lower than Esy^chiv) 
for either photon hehcity can therefore be excluded [63] . This is always a hard 
photon threshold, since the soft photon threshold occurs when the electron 
group velocity reaches the low energy speed of light, whereas the velocity 
required to produce any finite synchrotron frequency is smaller than this. 
The corresponding constraint is shown by the dashed line labelled "Synch. 
Cerenkov" in Fig. 1. This constraint improves on the current birefringence 




(23) 
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limit on ^ in the parameter region r] < 10^'^. At least one of the four pairs 
(±?7±, 1^1) must satisfy both this constraint and the synchrotron constraint 
discussed above. This amounts to saying that one of the two pairs {\r)±\, \^\) 
must satisfy this combined constraint or its reflection about the ^ axis. As with 
the hard Cerenkov constraint discussed above, this imposes a lower bound on 
one of |77±| once |^| is large enough. 

The synchrotron and Cerenkov constraints can further be linked, to obtain an 
upper bound on one of |r7±| as follows. We know that (automatically) at least 
one of the four parameters ±r)± satisfies the synchrotron constraint (a nega- 
tive lower bound), and at least one of them (automatically) satisfies the IC 
Cerenkov constraint (a positive upper bound). In fact, at least one of the four 
parameters must satisfy both of these constraints. The reason is that otherwise 
the synchrotron charges would violate the Cerenkov constraint, hence their 
energy would necessarily be under 50 TeV, which is thirty times lower than 
the Lorentz invariant value of 1500 TeV for the highest energy synchrotron 
charges. By itself this is not impossible since, as explained in Appendix A, 
if T] is positive a charge with lower energy can produce high frequency syn- 
chrotron radiation. However, the Crab spectrum is well accounted for with 
a single population of charges responsible for both the synchrotron radiation 
and the IC 7-rays. If there were enough charges to produce the observed syn- 
chrotron fiux with thirty times less energy per electron, then the charges that 
do satisfy the Cerenkov constraint would presumably be at least equally nu- 
merous (since their 77 is smaller so helicity decay and source effects would, 
if anything, produce more of them), and so would produce too many IC 7- 
rays [63]. Thus at least one of ±r]± must satisfy together the synchrotron, 
synchrotron-Cerenkov, and the IC Cerenkov constraints. That is, one of the 
four pairs {±r)±, ^) must fall within the dashed bell-shaped boundary in Fig. 1. 
This amounts to the statement that one of the pairs {r]±,C) must fall within 
the union of the bell-shaped boundary and its reflection about the ^ axis. This 
imposes the two-sided upper bound |77±| < 10"^ on one of the two parameters 
V±- 

5.9 Helicity decay 

The constraint Ar] = \r]^ — r]-\ < 4 on the difference between the LV param- 
eters for the two electron helicities was deduced by Myers and Pospelov [45] 
using a previous spin-polarized torsion pendulum experiment [87]. ^ Using the 
photon decay bound \rj±\ < 0.2 discussed here, we infer the stronger bound 

^ They also determined a numerically stronger constraint using nuclear spins, how- 
ever this involves four different LV parameters, one for tlie photon, one for the 
up-down quark doublet, and one each for the right handed up and down quark 
singlets. It also requires a model of nuclear structure. 
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Ary < 0.4. In this section we discuss the possibihty of improving on this con- 
straint using the process of hehcity decay. 

If 77- > 77+, negative helicity electrons are unstable to decay into positive helic- 
ity electrons via photon emission. (We assume for this section that rj^ > V+j 
the opposite case works similarly. ) This reaction has no kinematic thresh- 
old. However, the rate is very small at energies below an effective threshold 
(m^M/Ary)^/^ ^ 10 TeV (see Appendix C.7 for explicit rate calculations). The 
decay lifetime is minimized at the effective threshold. Below that it is longer 
than at least ~ (A?7)~^(p/10TeV)~^ x 10~^ seconds, while above it is given 
by ~ (p/lOTeV) x 10~^ seconds, independent of A77. 

While accelerator energies are well below the effective threshold if Arj is 0(1), 
in principle one might still get a bound by looking at fractional loss for a large 
population of polarized accelerator electrons. In practice this is probably im- 
possible due to other polarizing and depolarizing effects in storage rings. But 
just to see what would be required to improve the current bound suppose (op- 
timistically) that flipping of one percent of the electrons stored for 10"^ seconds 
could in be detected. The experimental exclusion of this phenomenon would 
require that the hfetime be greater than 10^ seconds. Using our overestimate of 
the decay rate (C.30) this would yield the constraint Arj < 10^(p/10 GeV)~^/^. 
To improve on the photon decay constraint would then require electron en- 
ergies at least around 200 GeV. When it was running LEPIl produced 100 
GeV beams, but the currently highest energy electrons in storage rings are 
around 30 GeV at HERA and 10 GeV at BABAR and the KEK B factory. 
Hence even if the non-LV polarization effects could be somehow factored out, a 
useful helicity decay bound from accelerators does not appear to be currently 
attainable. 

A helicity decay bound can also be inferred using the charged leptons in the 
Crab nebula. We have just argued that at least one of the two parameters 
?7± must lie within the dashed, bell-shaped region together with its reflection 
about the ^ axis. Let us assume that the pair [f]-,^) is inside the allowed 
region. We can divide the region into three parts, A where > 7 ■ 10^^, B 
where |r7_| < 7 • 10~^, and C where rj^ < —7 • 10~^. In the first case negative 
helicity electrons can produce the synchrotron radiation, but positive helicity 
positrons (which have LV parameter — ?7-) are below the synchrotron bound, 
so cannot. Now, if 77+ < —0.01 neither positive helicity electrons or negative 
helicity positrons can be responsible for the observed synchrotron radiation, 
since the electrons cannot emit the necessary frequencies and the positrons 
lose energy too rapidly via vacuum Cerenkov emission. Thus, in this case, of 
the four populations of leptons only the negative helicity electrons can produce 
the synchrotron radiation. However, an rj^ this low would allow rapid helicity 
decay (since the particle energy must be above 50 TeV which is around the 
effective hehcity decay threshold when A77 = 0.01) from negative to positive 
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helicity electrons, leaving no charges to produce the synchrotron. So if rj- 
is in region A (which violates Lorentz symmetry) we infer the lower bound 
?7_|_ > —0.01. Similarly, if r)_ is in region C, the same logic implies the upper 
bound 77+ < 0.01. If rj- is in region B, then negative hehcity electrons and 
positive helicity positrons are both able to produce the synchrotron radiation. 
Whatever the value of 77+, at least one of ±(?7- — 1]^) is positive, hence at 
least one of these two species is always stable to helicity decay, so no helicity 
bound can be presently inferred with one parameter in region B. If the Crab 
spectrum could be modelled and observed precisely enough to know that both 
(or all four) species must contribute to the synchrotron radiation, a helicity 
bound could be obtained in case B. 



6 Other types of constraints on LV 

In this section, we briefly summarize and point to references for various con- 
straints on LV effects besides those associated with 0{E/M) effects in QED. 



6. 1 Constraints on dimension 3 and 4 operators 

For the n — 2 term in (2,3), the absence of a strong threshold effect yields 
a constraint 772 ^ {m/pY{M/ n). If we consider protons and put ji = m = 
rUp ~ 1 GeV, this gives an order unity constraint when p ~ \fmM ~ 10^^ 
eV. Thus the GZK threshold (see the following subsection), if confirmed, can 
give an order unity constraint, but multi-TeV astrophysics yields much weaker 
constraints. The strongest laboratory constraints on dimension three and four 
operators for fermions come from clock comparison experiments using noble 
gas masers [88]. The constraints limit a combination of the coefficients for 
dimension three and four operators for the neutron to be below 10~^^ GeV 
(the dimension four coefficients are weighted by the neutron mass, yielding a 
constraint in units of energy). This corresponds to a bound on r]\ of order 10~^^ 
in the parametrization of (3) with \i = \ GeV. For more on such constraints 
see e.g. [33,89]. Astrophysical limits on photon vacuum birefringence give a 
bound on the coefficients of dimension four operators of 10~^^ [80]. 

6.2 Constraints at 0{E/M) from UHE cosmic rays 

In collisions of ultra high energy (UHE) protons with cosmic microwave back- 
ground (CMB) photons there can be sufficient energy in the center of mass 
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frame to create a pion, leading to the reaction p + 7cmb — ^ p + 7r°. The so 
caUed GZK threshold [34] for the proton energy in this process is 

Egzk ^ ^ ^ 3 X 10^°eV x (^^) (24) 



To get a definite number we have put equal to the energy of a photon at the 
CMB temperature, 2.7K, but of course the CMB contains photons of higher 
energy. This process degrades the initial proton energy with an attenuation 
length of about 50 Mpc. Since plausible astrophysical sources for UHE particles 
are located at distances larger than 50 Mpc, one expects a cutoff in the cosmic 
ray proton energy spectrum, which occurs at around 5 x 10^^ eV, depending 
on the distribution of sources [90] . 

One of the experiments measuring the UHE cosmic ray spectrum, the AGASA 
experiment, has not seen the cutoff. An analysis [91] from January 2003 con- 
cluded that the cutoff was absent at the 2.5 sigma level, while another experi- 
ment, HiRes, is consistent with the cutoff but at a lower confidence level. (For a 
brief review of the data see [90].) The question should be answered in the near 
future by the AUGER observatory, a combined array of 1600 water Cerenkov 
detectors and 24 telescopic air fluorescence detectors under construction on 
the Argentine pampas [92]. The new observatory will see an event rate one 
hundred times higher, with better systematics. 

Many ideas have been put forward to explain the possible absence of the 
GZK cutoff [90], one being Lorentz violation. According to equation (24) the 
Lorentz invariant threshold is proportional to the proton mass. Thus any LV 
term added to the proton dispersion relation E"^ = p'^ + will modify the 
threshold if it is comparable to or greater than rii^ at around the energy 
Egzk- Modifying the proton and pion dispersion relations, the threshold can 
be lowered, raised, or removed entirely, or even an upper threshold where the 
reaction cuts off could be introduced (see e.g. [53] and references therein). 

If ultra-high energy cosmic rays (UHECR) are (as commonly assumed) pro- 
tons, then strong constraints on n = 3 type dispersion can be deduced from a) 
the absence of a vacuum Cerenkov effect at GZK energies and b) the position 
of the GZK cutoff if it will be actually found. 

a) For a soft emitted photon with a long wavelength, the partonic structure 
of a UHECR proton is presumably irrelevant. In this case we can treat the 
proton as a point particle as in the QED analysis. With a GZK proton of 
energy 5 x 10^^ eV the constraint from the absence of a vacuum Cerenkov 
effect is ?7 < 0{\Q~^*) [53]. Since the helicity of cosmic rays is not observed, 
one can say only that this constraint must be satisfled for at least one helicity. 
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For a hard emitted photon, the partonic nature of the proton is important 
and the relevant mass scale will involve the quark mass. The exact calculation 
considering the partonic structure for n = 3 has not been performed. Treating 
the proton as structureless, the threshold region would be similar to that 
in [53] . The allowed region in the rj — ^ plane would be bounded on the right 
by the ^ axis (within a few orders of magnitude of 10~^^) and below by the 
line ^ = 1] [53]. This constraint applies to both photon helicities, but only to 
one proton helicity, since the UHECR could consist all of a single helicity. In 
principle, however, what one can really constrain is some combinations of the 
various quark dispersion parameters. This approach has been worked out in 
detail using parton distribution functions in [93]. 

b) If the GZK cutoff is observed in its predicted place, this will place limits on 
the proton and pion parameters rjp and r/^. For example, if the GZK cutoff is 
eventually observed to be somewhere between 2 and 7 times 10^^ eV then there 
are strong constraints of O(10~^^) on the relevant rjp and r^jr [53]. (Allowing for 
helicity dependence, no set of parameters allowing long distance propagation 
(forbidding the vacuum Cerenkov effect) should modify the GZK cutoff.) 

As a final comment, an interesting possible consequence of LV is that with 
upper thresholds, one could possibly reconcile the AGASA and Hi-Res/Fly's 
Eye experiments. Namely, one can place an upper threshold below 10^^ eV 
while keeping the GZK threshold near 5 x 10^*^ eV. Then the cutoff would 
be "seen" at lower energies but extra flux would still be present at energies 
above 10^*^ eV, potentially explaining the AGASA results [53]. The region of 
parameter space for this scenario is terribly small, however, again of O(10~^^). 



6.3 Constraints on dimension 6 operators 

As previously mentioned, CPT symmetry alone could exclude the dimension 
five LV operators in QED that give 0{E/M) modifications to particle disper- 
sion relations. Moreover, the constraints on those have become quite strong. 
Hence we close with a brief discussion of the constraints that might be possible 
at 0{E'^ / M"^). Such LV effects arise for example from dimension 6 operators. 
Note that helicity dependence of the LV parameters is not required in this 
case, and on the other hand it can occur without violating CPT. Without any 
particular theoretical prejudice, one should keep in mind that constraints will 
generally only limit the parameters for some helicity species, while they might 
be evaded for other hehcities. 

As previously discussed we can estimate that the LV modification of the dis- 
persion relation becomes important when comparable with the mass term. 
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r]4p'^/M'^ < w?, which yields 




m lO^^eV 
IMeV p 
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(25) 



Thus, for electrons, an energy around 10^" eV is needed for an order unity 
constraint on 774, and we are probably not going to see any effects directly 
from such electrons. 

For protons an energy ~ 10^^ eV is needed. This is well below the UHE cosmic 
ray energy cutoff, hence if and when Auger [92] confirms the identity of UHE 
cosmic rays as protons at the GZK cutoff, an impressive constraint of order 
774 ^ 10~^ will follow from the absence of vacuum Cerenkov radiation for 10^'^ 
eV protons. From the fact that the GZK threshold is not shifted, a constraint 
of order 774 > — 10~^ will follow, assuming equal 774 values for proton and pion. 

In fact, if one assumes the cosmic rays already observed near but below the 
GZK cutoff are hadrons, one obtains a strong bound [93]. Depending on the 
species and helicity dependence of the LV coefficients, bounds of order 10~^ or 
better can be placed on 774. The bounds claimed in [93] are actually two sided, 
and come about in a manner somewhat analogous to (but more complicated 
than) the photon decay constraint (see section 4.3.1). They are derived by 
using a parton model for particles where the LV coefficients apply to the 
constituent partons. By considering many different outgoing particle spectra 
from an incoming hadron in combination with the parton approach the authors 
of [93] are able to find sets of reactions that yield two sided bounds. Hence, 
the parton approach is quite useful, as it dramatically increases the number 
of constraints that can be derived from a single incoming particle. 

One might think that impressive constraints can also be obtained from the 
absence of neutrino vacuum Cerenkov radiation: putting in 1 cV for the mass 
in (25) yields an order unity constraint from 100 TeV neutrinos, but only if the 
Cerenkov rate is high enough. The rate will be low, since it proceeds only via 
the non-local charge structure of the neutrino. Recent calculations [94] have 
shown that the rate is not high enough at that energy, even for cosmogenic 
neutrinos. However, for lO^'' eV UHE neutrinos, which may be observed by 
the proposed EUSO (Extreme Universe Space Observatory) [95] and/or OWL 
(Orbiting Wide-angle Light collectors) [96] satellite observatories, the rate 
will be high enough to derive a strong constraint as long as the neutrinos are 
cosmogenic, and perhaps even if they originate closer to the earth. The specific 
value of the constraint would depend on the exact value of the rate, which 
has not yet been computed. For a gravitational Cerenkov reaction, the rate 
(which is lower but easier to compute than the electromagnetic rate) would be 
high enough for a 10^° eV neutrino from a distant source to radiate provided 
7^4 > 10~^. Hence in this case one might obtain a constraint of order 7^4 < 10~^ 
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from gravitational Ccrcnkov. If Ccrenkov constraints apply only when the 
observed UHE neutrino originates from a distant source, one would need to 
either identify an astrophysical object as the source or somehow otherwise rule 
out local generation of the neutrino. 

A time of flight constraint at order (E/M)'^ might be possible [97] if gamma 
ray bursts produce UHE (~ 10^^ eV) neutrinos, as some models predict, via 
limits on time of arrival differences of such UHE neutrinos vs. soft photons (or 
gravitational waves). Another possibility is to obtain a vacuum birefringence 
constraint with higher energy photons [81] , although such a constraint would 
be less powerful since EFT does not imply that the parameters for opposite 
polarizations are opposite at order (E/M)"^. If future GRB's are found to 
be polarized at ~ 100 MeV, that could provide a birefringence constraint 

|e4+-C4-| <1. 



7 Future prospects 

In the last decade or so the old dogma that all our observations are insensitive 
to Planck scale effects has been shown to be quite wrong if Lorentz symmetry is 
violated. A large number of impressively strong constraints have been obtained 
on LV Planck scale effects arising from operators of mass dimension three, four, 
five, and even six. 

Prom the conceptual point of view, the most burning issue is naturalness of 
small low energy LV. We have explained in this article both the rationale for 
adopting an effective field theory paramctrization of LV, and the "naturalness" 
problem of preserving approximate low energy Lorentz invariance if LV is 
to exist in the UV theory. As discussed in section 3.4, currently the best 
prospect for resolving this issue in favor of LV is via supersymmetry, but 
perhaps other ideas could work. In the Lorentz violating, softly broken SUSY 
QED framework of Ref. [61], dimension five (CPT violating) LV operators are 
already strongly constrained by current QED experiments, but the dimension 
six LV operators are not yet constrained at the Mpf^^^^. level. However the 
authors of that reference indicate that Planck level constraints on dimension 
six LV quark operators in the standard model may be possible. We emphasize 
that high energy LV dispersion of the sort discussed in the present paper does 
not occur in the scenario of Ref. [61]. 

Constraints continue to improve with advances in observational technology. In 

this concluding section wc examine some of the future prospects for improving 
the constraints on high energy LV using astrophysical observations. To begin, 
we note that while photon "time of flight" constraints are not competitive with 
other current constraints, they enjoy a special status in being less dependent 
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on assumptions aboTit the underlying theoretical framework. Given that time 
of flight constraints are strongest with a high energy signal that has structure 
on short time scales, gamma ray bursts are probably the best objects for im- 
proving these constraints. Experiments hke GLAST (Gamma-ray Large Area 
Space Telescope) [98] may be helpful in improving bounds from time of flight; 
they might lead to a better understanding of GRB emissions and provide far 
better photon counts and time resolution (less than lOfisec). These improve- 
ments might help overcome some of the problems involving GRB time of flight 
measurements discussed in section 5.1. See e.g. [78] for a detailed discussion 
of these issues. 

Here we focused on LV in the dimension five, CPT violating sector of QED, 
characterized by the photon parameter \^\ bounded at O(10~^) and the elec- 
tron parameters \r]±\ bounded collectively at O(10~^), while at least one of 
them is bounded at O{10~^). The most promising prospect for significant im- 
provement in these bounds would be a closer study of the effect of this sort of 
LV on the complete electromagnetic spectrum of the Crab nebula and or other 
plerions. We have been very conservative in imposing constraints, by allowing 
for the possibility that with LV the synchrotron radiation could be produced 
by much lower energy charges than required in the usual Lorentz invariant 
case. We have also tolerated the possibility of four different populations of 
charges behaving differently (two helicities for electrons and positrons). By 
modelling more completely the production of the electron-positron wind and 
electromagnetic spectrum in the presence of such LV it may be possible to 
infer that charges of energy 1500 TeV are in fact required as in the Lorentz 
invariant case. If so, the constraints on ri± might be improved by a factor of 
order ~ (50/1500)'^ ~ 10~"^ to \ri±\ < 10"^, and that on |,^| by a factor of order 
~ 10~^ to |,^| < 10^® as well. This would seem to be rather definitive so is a 
useful goal. To achieve it would perhaps require not only better modelling of 
the LV effects but more precise observations of the Crab nebula. In this re- 
gard we can look forward to the observations to be made with GLAST which 
should start its observations in 2007. 

The birefringence bound, currently of order |^| < 10^'' from UV polarization 
observations of a distant galaxy, might be improved in various ways. Since the 
LV effect scales as the square of the energy, the most dramatic improvement 
would be if polarized gamma rays from a distant source could be observed, 
for example from a GRB. (The observation reported in [82] appears now to 
be unreliable.) This might be possible with RHESSI or via the IBIS camera 
of the INTEGRAL satellite [99], using the so-called Compton mode that was 
used in the previous RHESSI analyses. Less dramatically, it might instead be 
possible to improve by an order of magnitude the constraint reported in [42], 
obtained using UV light from the radio galaxy 3C 256, via specifically targeted 
observations possibly using FORSl [100] at the ESO Very Large Telescope 
(VLT), and/or by observing more distant sources of polarized UV light. 
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Bounds on LV effects on threshold reactions for ultra-high energy cosmic rays 
can be significant even when suppressed by two powers of the Planck mass, i.e. 
at order E"^ /M"^. This is important since CPT symmetry alone could preclude 
effects at order E/M. We already discussed the currently available bounds. 
The possibility of improvement lies in experiments which will be capable of 
reaching high sensitivities at the small fluxes associated with the highest en- 
ergy particles. The Auger detector will be the next big thing, with first release 
of data expected in summer 2005. 

It is conceivable that ultra-high energy neutrinos could provide extremely 

stringent constraints on LV due to their small mass. The vacuum Ccrenkov 
effect occurs at a given energy with a minimum LV coefficient that scales 
with mass as rn?. On the other hand the rate for the reaction is very low, so 
it is not clear whether in practice this possibility will pan out. UHE neutri- 
nos might be observed using IceCube [101], or planned detectors such as the 
cited satelhtes EUSO [95] and OWL [96] (the former being already in an ad- 
vanced stage of development). Also promising are upcoming experiments like 
the sea based detector NESTOR (Neutrino Extended Submarine Telescope 
with Oceanographic Research) [102] which is now starting its activity, or the 
balloon-borne detector ANITA (Antarctic Impulse Transient Antenna) [103] 
whose first flight is planned for the end of 2006. Finally also very effective 
could be some proposed experiments like the underground salt based detector 
SalSA (Salt-dome Shower Array) [104]. 

In closing, we have been motivated by considerations of quantum gravity to 
look for evidence of, and bounds on, Lorentz violation. These considerations 
include suggestions of LV from various viewpoints, as well as simple despera- 
tion for observational guidance. The constraints discussed here and elsewhere 
are mounting up to a signiflcant limitation on the possibility of LV in quantum 
gravity, at least insofar as would appear in low energy effective fleld theory. 
It seems very worthwhile to push on, even if the ultimate result will only be 
extremely restrictive constraints. In this way we will have solidifled the obser- 
vational foundations of Lorentz symmetry, and acquired valuable constraints 
on future speculations on the nature of quantum gravity. 



A LV cutoff of synchrotron radiation frequency 



Cycling electrons in a magnetic fleld B emit synchrotron radiation with a 
spectrum that sharply cuts off at a frequency ujc given in the Lorentz invariant 
case by the formula 

^. = \eB'-^, (A.1) 
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where "y{E) = (1 — v^{E)/c^)~^/'^ . Here v{E) is the electron group velocity, 
and c is the usual low energy speed of light. In standard relativistic physics, 
E = 7m, so the energy dependence in (A.l) is entirely through the factor 7^, 
which grows without bound as the energy grows. Arbitrarily large synchrotron 
frequencies are therefore possible. 

In the Lorentz violating case with negative values of 1], electrons have a maxi- 
mal group velocity strictly less than the low energy speed of light, hence there 
is a maximal synchrotron frequency that can be produced, regardless of the 
electron energy [46].^° As we shall argue below, Eq. (A.l) still holds in this 
case, assuming the framework of effective field theory. Therefore the maximal 
frequency is obtained by maximizing '~f^{E)/E with respect to the electron 
energy. Using the difference of group velocities 



obtained from the electron dispersion relation (13), we find that this maxi- 
mization yields 

^max ^ Q_34 (A.3) 



This maximum frequency is attained at the energy i?max = {—2m'^M/5r]y^^ = 
10 {—rj)~^^^ TcV. This is higher than the energy that produces the same cutoff 
frequency in the Lorentz invariant case, but only by a factor of order unity. 

Note that if r) is positive, then the effect is the opposite: an electron can 
produce a given frequency of synchrotron radiation with an energy less than 
in the Lorentz invariant case. The electron speed can even exceed the low 
frequency speed of light, at which point ^{E) diverges. This corresponds to 
the soft Cerenkov threshold discussed in section C.5. 

We shall now justify the derivation of (A.l) in the LV case, adding some details 
to previously published work [46]. The first step is a purely kinematical analysis 
that does not assume Lorentz invariance, which follows the standard heuristic 
derivation [105]. A cycling electron of energy E emits radiation in a cone of 
some opening angle S{E). The cone sweeps past a distant observer as the 
electron moves on a circle of radius R{E) through an angle S{E). The electron 
travels at a speed (group velocity) v{E) so the time it takes to orbit through 

Since we have LV in the photon sector as well, there is no "speed of light" per 
se. However, the emitted frequencies of synchrotron radiation are much lower than 
the energy of the source particles in the Lorentz invariant case. Effectively the LV 
violation in the photon sector can be ignored. 
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the angle S{E) is At = R{E)6{E)/v{E). The hght from the leading edge 
of the cone travels a distance c{uj)At while the electron travels the distance 
v{E)At toward the observer. Hence the spatial width of the pulse seen by the 
observer is approximately {c{uj) — v{E j)At, which arrives at the observer over 
a time interval equal to this distance divided by the speed of light. The cut off 
frequency of the synchrotron pulse is roughly the inverse of this time interval, 

^ 3 1 1 

AR{E)S{E) c{uj,)-v{E)' ^ 



where we have used the fact that the electron and photon speeds are very close 
to the low energy speed of light c, which is set equal to unity. In the Lorentz 
invariant case the radius is given by R{E) = E/eB and the opening angle is 
S{E) ~ 7~^(£'). The numerical constant in (A. 4) is chosen so that when these 
values are substituted the correct relativistic result (A.l) is obtained. 

Under the assumption that the Lorentz violation is described by effective field 
theory, we will argue that R{E) and 5{E) are the same in the LV case as in 
the Lorentz invariant case. Moreover, since the emitted photons have relatively 
low energy compared to the electrons, it turns out that ^ can be neglected in 
the relevant region of parameter space. (As shown in [46] the photon energy is 
low enough to neglect any possible LV correction as long as |^| < 10-'^^(— 77)^/^.) 
Thus c{uJc) in (A. 4) can be replaced by c, so the reciprocal of the difference of 
group velocities is well approximated by 2^'^{E), as in the Lorentz invariant 
case. This yields (A.l) where now all the dependence on the LV is in how the 
gamma factor depends on E. 

The radius R{E) is determined by the equation of motion for the electron 
in a magnetic field. All Lorentz violating terms in the equation of motion 
are suppressed by the ratio E/M. From this it is clear that the trajectory 
determined by a given initial position and momentum is not much affected 
by the Lorentz violation, as long as E remains much smaller than the Planck 
energy. 

This conclusion about R{E) can be more explicitly verified by reference to 
the modified Hamilton's equations, with the Hamiltonian given in terms of 
the momentum by the dispersion relation between energy and momentum. 
The leading high energy corrections come from modifications to the minimal 
coupling terms. As usual the minimal coupling is incorporated replacing the 
momentum by p — eA, where A is a vector potential for the magnetic field. 



A subtlety here is that this argument fails if wc refer to the velocity rather than 
to the momentum. The reason is that tiny differences in the speed make a huge 
difference in the momentum and energy when the speed is close to the speed of 
light. 
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This yields the equation of motion a = [1 + 3r]E/2M](e/ E) v x B, where 
we have kept only the lowest order term in 77 and assumed relativistic energy 
E ^ m. Since E <^ M, the presence of the Lorentz violation makes very 
little difference to the orbital equation, hence we conclude that to a very good 
approximation the radius is related to the magnetic field and the energy of 
the electron by the standard formula R{E) = E/eB (where again the speed 
of the electron has been set equal to unity). This result was independently 
confirmed by [106] from the starting point of the Dirac equation. 

The angle S{E) scales in the Lorentz invariant case as ^~^{E). Since R{E) and 
hence the charge current is nearly unaffected by the Lorentz violation, any 
significant deviation in S{E) could only come from the modified response of 
the electromagnetic field to the given current. The Lagrangian in the presence 
of the dimension five LV operator for the EM field in (6) takes the form 

L = ^F^.F^" + j^G(A) + AJ% (A.5) 



where ^G{A)/M is the LV operator and is the given current. The equations 
of motion are 



If we expand A" = + ^i, where is the field in the Lorentz invariant 
case, then to lowest order the LV correction A^ satisfies 



The correction A^ is suppressed by at least one power of M. In terms of Fourier 
components, this suggests that the relative size of the LV correction is of order 
^Lv/M, which is extremely small for even the peak photon energies ~ 10^ eV 
present in the synchrotron radiation from the Crab nebula. Therefore near 
the source current the amplitude of the electromagnetic field modification is 
negligible in comparison to the zeroth order synchrotron radiation (which is 
still present), and hence S{E) scales with ^^^{E) in the usual way. Note that 
this does not mean that the relative size of the correction is small everywhere. 
The source of A^ on the right hand side of (A. 7) is non- vanishing not just 
where the current lies but also wherever the emitted radiation propagates. 
It can therefore produce a cumulative effect that can be large. In fact, this 
must clearly happen, since the group velocity is modified by the LV term, so 
that after enough propagation the relative correction must become of order 
unity so that it can shift the support of the radiation to a disjoint location. 
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A detailed analysis of synchrotron radiation in this LV QED theory was carried 
out in Ref. [107]. Our heuristically obtained results were confirmed there, and 
certain regimes were identified in which the LV solution deviates strongly from 
the Lorentz invariant one. Our interpretation of such deviations is what was 
mentioned at the end of the previous paragraph, i.e. that these are cumulative 
effects. A claim of significant deviations was made in Ref. [108] however we 
do not currently understand the basis of that claim. 



B Threshold configurations 

Threshold configurations and new phenomena in the presence of LV disper- 
sion relations were systematically investigated in [39,69] (see also references 
therein). We give here a brief summary of the results. We shall consider reac- 
tions with two initial and two final particles (results for reactions with only one 
incoming or outgoing particle can be obtained as special cases) . Following our 
previous choice of EFT we allow each particle to have an independent disper- 
sion relation of the form (1) with E{p) a monotonically increasing non- negative 
function of the magnitude p of the 3-momentum p. While the assumption of 
monotonicity could perhaps be violated at the Planck scale, it is satisfied for 
any reasonable low energy expansion of a LV dispersion relation. EFT further 
implies that energy and momentum are additive for multiple particles, and 
conserved. 

Consider a four-particle interaction where a target particle of 3-momentum 
P2 is hit by a particle of 3-momentum pi, with an angle a between the two 
momenta, producing two particles of momenta ps and p4. We call j3 the angle 
between ps and the total incoming 3-momentum pj„ = pi + P2. We define the 
notion of a threshold relative to a fixed value of the magnitude of the target 
momentum p2- A lower or upper threshold corresponds to a value of pi (or 
equivalently the energy Ei) above which the reaction starts or stops being 
allowed by energy and momentum conservation. 

We now introduce a graphical interpretation of the energy-momentum con- 
servation equation that allows the properties of thresholds to be easily under- 
stood. For given values of (pi,p2,a, /3,P3), momentum conservation determines 
P4. Since p^ and p^ determine the final energies E3 and i?4, we can thus define 
the final energy function E'^'^'^'^(j)i). (Since p2 is fixed we drop it from the 
labelling.) Energy conservation requires that Ef be equal to Ei{pi), the initial 
energy (again, we do not indicate the dependence on the fixed momentum P2). 

Now consider the region TZ in the {E,pi) plane covered by plotting i?J'^'^^(pi) 
for all possible configurations {a,(3,p3). An example is shown in Figure B.l. 
The region TZ is bounded below by £■ = since the particle energies are 
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Fig. B.l. Graphical representation of energy- momentum conservation in a 
two-particle reaction. TZ is the region covered by all final energy curves E'^'^'^^{pi) 
for some fixed p2, assuming momentum conservation holds to determine ^4. The 
curve Ei{pi) is the initial energy for the same fixed p2- Where the latter curve lies 
inside TZ there is a solution to the energy and momentum conservation equations. 
In the example shown there is both a lower and an upper threshold for the reaction. 

assumed non- negative, hence it has some bounding curve Eb{pi). Similarly 
one can plot Ei{pi). The reaction threshold occurs when Ei{pi) enters or 
leaves TZ, since it is precisely in TZ that there is a solution to the energy and 
momentum conservation equations. 

This graphical representation demonstrates that in any threshold configura- 
tion (lower or upper) occurring at some pi, the parameters {a,f3,p^) are such 
that the final energy function E'^'^'^'^{pi) is minimized. That is, the configura- 
tion always yields the minimum final particle energy configuration conserving 
momentum at fixed pi and p2- Prom this fact, it is easy to deduce two general 
properties of these configurations: 

(1) All thresholds for processes with two outgoing particles occur at parallel 
final momenta {(3 = 0). 

(2) For a two particle initial state the momenta are antiparallel at threshold 
{a — tt). 

These properties are in agreement with the well known case of Lorentz in- 
variant kinematics. Nevertheless, LV thresholds can exhibit new features not 
present in the Lorentz invariant theory, in particular upper thresholds, and 
asymmetric pair creation. 
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Figure B.l clearly shows that LV allows for a reaction to not only to start 
at some lower threshold but also to end at some upper threshold where the 
curve Ei exits the region TZ. It can even happen that enters and exits TZ 
more than once, in which case there are what one might call "local" lower and 
upper thresholds. 

Another interesting novelty is the possibility to have a (lower or upper) thresh- 
old for pair creation with an unequal partition of the initial momentum 
into the two outgoing particles (i.e. Ps ^ P4 ^ pin/2). Equal partition of mo- 
mentum is a familiar result of Lorentz invariant physics, which follows from 
the fact that the final particles are all at rest in the zero-momentum frame at 
threshold. This has often been (erroneously) presumed to hold as well in the 
presence of LV dispersion relations. 

A reason for the occurrence of asymmetric LV thresholds can be seen graph- 
ically, as shown in Figure B.2. Suppose the dispersion relation for a massive 
outgoing particle E^utip) has negative curvature at p = Pinf^, as might be 
the case for negative LV coefficients. Then a small momentum-conserving dis- 
placement from a symmetric configuration can lead to a net decrease in the 
final state energy. According to the result established above, the symmetric 
configuration cannot be the threshold one in such A lower pi could 

satisfy both energy and momentum conservation with an asymmetric final 
configuration. A sufficient condition for the pair-creation threshold configura- 



AE,(p,) 




P3= P4=Pln/2 



Fig. B.2. Asymmetric pair production. The negative curvature of the outgoing 
particle dispersion relation allows to save energy by providing the pair partners 
with different portions of the initial momentum 

tion to be asymmetric is that the final particle dispersion relation has negative 
curvature at p = Pm/2. This condition is not necessary however, since it could 
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happen that the energy is locally but not globally minimized by the symmetric 
configuration. 



C Rates and thresholds of LV processes 



In this appendix we first derive a general expression for reaction rates. We 
then discuss the rates and thresholds for three reactions of interest: photon 
decay, vacuum Cerenkov, and helicity decay. 



C.l General rate expressions 



The first step in deriving constraints from particle interactions is to derive 
the rate F for various processes. As it turns out, the most useful constraints 
are derived from reactions where a single particle decays into two particles. 
If we label the incoming particle momentum and helicity/polarization by p, s 
and the final particles by p', s', p", s" then the rate for fixed helicities is given 
by [109] 



nP.s.s\s")^ [^ .,JZ%^ |M(p,.,p-,.-,p-,.-)r (C.l) 

xd{Ep,, - e;, - E;.^,,)d^'\p - p' - p") 

where M(p, s, p', s', p", s") is the matrix element. As mentioned in section 
4.1.2 our fermion spinors are not quite normalized to 2E. As a result (C.l) 
is not quite correct. However, as long as we work in the preferred frame, this 
correction is small and can be neglected since it is only a factor in overall 
normalization. 

If we choose the initial particle to be travelling in the ^-direction, then we can 
integrate over p" and apply the axial symmetry to get 

r(p, .,.',.") = / -^ j^ y^^^ \MfsiEp, - E'^,, - e;_^,mc.2) 

J l&TT Ep^sEp,^,,Ep_p,^^„ 



Integrating over then yields 



IMI 



p',g" 



. (C.3) 
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This expression and those that follow are evaluated at the solution of the 
energy conservation equation, p'j_ = p'±{p'z,p, s, s', s"). The limits of integration 
P'ziiP'z2 the bounds between which such a solution exists. We are mostly 
interested in cases where there is no solution unless p is greater than some 
threshold momentum ptu- 

For reactions that do not occur in the Lorentz-invariant limit, the phase space 
must close up as we approach Lorentz invariance, i.e. the magnitude of p'^ 
vanishes as Epi — > oo. Hence p'^ must be small relative to p'^ and p — p'^ if 
p <^ Epi. The exception to this is if either p^jP — p'^ is near zero, i.e. when 
either or p'.^2 ~ P- This region is, however, a negligible amount of phase 

space and so we ignore it. Thus p'^ can be treated as an expansion parameter. 



Expanding the derivative of E'^t^^i + E'^-p'^s" ^^h respect to p'j^ gives 



-1 T?l 77" 



(C.4) 



p'±Ep^s 

and substituting into (C.3) gives 

r(p,.,.',.'')^/jijf |M|^. (C.5) 



zl 



To compute an energy loss rate for the Cerenkov effect a p'^ must be inserted 
in the integrand to reflect the energy carried off by the emitted photon. 



C.2 Photon decay matrix element 

We now show that, if the incoming momentum p is above the threshold pth, 
then the rate for photon decay is very high. 

The QED matrix element for photon decay is 

iM = ieus{p')e^i''vs{p") (C.6) 

where u{p') is the outgoing electron spinor and v{p") is the outgoing positron 
spinor. We consider the case where the electron and positron have positive and 
negative helicity respectively, so only the coefficient t]^ is involved. The LV 
term in the electron dispersion is r^+p^ (with m"*" wavefunction (12)) while that 
in the positron dispersion is —rj^p^ (with v~ wavefunction (16)). The result 
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is directly transferable to the rj- case. We will also relabel E{p) as u!{k) since 
the incoming particle is a photon. 



With these choices the dominant contribution to the matrix element at high 
energy is straightforwardly evaluated as 



iM = ieV2E^V2E^ [x+(p')He± ■ ^)x+ip") 



(C.7) 



which is equal to 
ieV2E'E" 




sin [^) (1±1) - sin 



cos 




(C.8) 



where 6' , 9" are the opening angles of the electron and positron and the ± 
reflects the initial photon polarization. (The overall phase depends on the 
phase convention for the spinors and does not affect the result.) Note that the 
matrix element vanishes in the threshold configuration when the momenta are 
all parallel. This is because wc chose the electron and positron helicities to be 
opposite. Nevertheless, above threshold the rate is high enough to obtain a 
useful constraint, as we now argue. 

The opening angles are small, since the perpendicular momenta are Planck 
suppressed. Therefore we can expand (C.8) to first order in 6*', 6*", yielding 

iM = \esjE'E"l2 



4(1 ±1)_ 4(1^1) 



(C.9) 



where j)'^ is the magnitude of the perpendicular momentum (which is the same 
for the electron and positron). Recall that "p'^ is determined from the energy 
conservation equation (where k is the initial photon momentum), 

M p'^ p'^ M k-p'^ k-p'^ M ■ ' 



In (CIO) we have already cancelled the 0th order terms using momentum 
conservation, and we have neglected higher order terms in p'^/pz- We have 
also assumed that the z-components of the electron and positron momenta 
are both positive. This is not necessarily the case, but when m <^ k <^ M 
(which has also been assumed) we are only missing a tiny fraction of the phase 
space by making this restriction. 

Defining z — p'^/k e [0, 1], p'^ is given by 

pI = [H - vA'^z - l)]z{l -z)^- m\ (C.ll) 
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If the reaction is to happen at all, then the RHS of (C.ll) must be positive. 
The points where = give the upper and lower bounds 'p'z\^v'z2 for the 
integral (C.5). Well above threshold the mass term becomes negligible, and 
the longitudinal phase space is determined by those z values for which ±^ — 
77+ (22; — 1) is positive in the interval z e [0, 1]. 

Using these results in the expression (C.5) for the decay rate we find that, well 
above threshold, the rate for decay of a positive helicity photon to a positive 
helicity electron and negative helicity positron is given by 



2iL2 



487rM 



(2^ + ^+) e>(e + ^+) + %^^e(|77+| 



(C.12) 



where ©(a;) is the Heaviside step function. This agrees up to a numerical factor 
with the estimate in [67] in the case where 77+ = 0. 

For the purposes of imposing constraints, it is important to know how rapidly 
the decay rate becomes large above threshold. We will return to this question 
in section C.4, after we derive the photon decay threshold. 



C.3 Photon decay threshold 



The photon decay threshold is found by setting p'l — and solving for the 
minimum value of k for which there is a solution for some z. Prom the thresh- 
old analysis described in section B we know that the outgoing momenta are 
parallel, hence < ;2 < 1 at threshold, so the minimum must have z within 
this range. 

Let us specialize to the case of a positive helicity photon (the negative helicity 
case works similarly). In terms of the variable x = 1 — 2z, the condition that 
(C.ll) vanish then becomes 

(l-x2)(C + r/+x) = — ^. (C.13) 



The range of x is —1 < x < 1, so we see that if \ri^\ is sufficiently large 
for a given ^ the conservation equation has a solution. Hence the threshold 
depends on 77+ only through its absolute value (the equation is invariant under 
r7+,x — > — 77+,— x). Physically, this corresponds to the fact that since the 
positron and electron have opposite dispersion, one is always subluminal. By 
depositing most of the outgoing momentum in the subluminal particle we can 
magnify its LV while minimizing the LV of the superluminal partner. Hence 
for either sign of 77+ we can find a solution. 



45 



Since the RHS of (C.13) decreases as k increases, the threshold will occur 
when the LHS is at its maximum in x. Let us consider a few special cases. 
First, if If r;_|_ = 0, then the LHS has a positive maximum only if ^ > 0. This 
occurs at x = (equal electron and positron momenta) and is equal to ^. In 
this case the threshold kth satisfies ^ = Arri^M/k^f^. If instead ^ = 0, then the 
maximum of the LHS occurs at a; = l/v^S, and is equal to 2|?7+|/3v^- In this 
case the threshold kff,, satisfies \r]^\ = G^/Srn^ M / kfj.^. If < — the LHS is 
never positive in the range — 1 < a; < 1, so photon decay does not occur for 
such parameters. More generally, the maximum of the LHS occurs at 

•^max o 



Substituting Xmax into (C.13) allows us to find the maximum of the LHS, but 
the expression is a bit complicated so we don't bother to display it here. One 
can use it to solve for ^ as a function of 77+ and kth and thus derive a constraint. 



C.4 Photon decay rate near threshold 

We now return to the question of reaction rate above threshold. If particle 
lifetimes are short slightly above threshold, then we are justified in using the 
absence of a threshold to establish constraints on the parameters. As a simple 
example of the rapidity with which reaction rates increase above threshold, 
let us analyze photon decay with ^ = 0, a positive helicity incoming photon, 
and 7^+ > (other parameter choices and the Cerenkov effect work similarly) . 

We expand the initial photon momentum as A; = (1 + a)kth, where kff^^ — 
6-\/3m^M/|r7+| and a <S 1. At threshold the entire phase space consists of 
one point, x — x^ax — I/a/S, and p'^ vanishes at this point. To calculate the 
rate using Eqns. C.5 and C.9 we need to know the range in x slightly above 
threshold and the value of p'^ in this range to first order in a. The rate is then 
given by 

Xmax ^X 

where wc have rewritten everything in terms of x. Ax is the spread around 
Xmax as wc uiovc abovc threshold. The fact that x^ax is an extremum implies 
that spread around x^ax is symmetric, and just above threshold the spread is 
small. Ax <C Xmax- Expanding pi about Xmax as a parabola, we find Eq. C.15 



(C.14) 



46 



is approximately given by 

m = '\ 1+^— p'j(w)Ax. (C.16) 



Now p^(x^ax) is given by Eq. C.ll at A; = kth{l + a), i.e. 

1 - 



From the definition of kth we know tliat tlic zcrotli order term in a vanislies. 
Witli tliis fact we can easily solve for pl{xmax) — "ioLm^ (to first order in a). 

All that remains to evaluate is Ax. At the endpoints of integration in Eq. C.5 
we know that = 0, which implies that we have the relation 

r7+(x„,,. ± Aa;)i^^^^^^^^^^A;,\(l + af - = 0. (C.18) 



Expanding Eq. C.18 to lowest order in Ax and a yields 

2 2 2 4:771^ 

Xmaxi^ ~ X^^^) + SaXjnaxi^ ~ X^g^^) — Sx^axi^x) ~ ^ ■ (C.19) 



There is no linear term in Ax since x^ax is an extremum of Eq. C.18. The 

first term on the left hand side is equal to the right hand side by definition of 
kth so we obtain (Ax)^ = 2a/3. Substituting Ax and pl{xmax) into Eq. C.16 
yields our final expression for the rate as a function of a, 

r = p^J^Z(2 + x/3)a3/2 ^C.20) 



To see how quickly the reaction starts to happen above threshold, consider 
an incoming photon of energy 1% above threshold. The lifetime of the photon 
is equal to 10~^/ct/i/10 TeV seconds, short enough that on any relevant astro- 
physical timescale a population of photons of multi-TeV energies will almost 
completely decay. As this example demonstrates, since the lifetime is extremely 
short only slightly above threshold, we are justified in using threshold values 
to derive constraints. 
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C.5 Vacuum Cerenkov threshold and rate 



The vacuum Cerenkov process e — > 67 is forbidden by angular momentum 
conservation in the threshold configuration. However, as with photon decay 
to opposite hehcities, the rate becomes large above threshold. In fact the 
rate calculation is very close to that for photon decay. The only significant 
difference is that a factor of the outgoing photon energy uj{k) must be inserted 
into (C.5) since wc are concerned with the energy loss rate. The net result is 
that dE/dt ~ p^/M above threshold, which imphes that a 10 TeV electron 
would emit a significant fraction of its energy in 10~^ seconds. Hence in this 
case as well a threshold analysis will give accurate results. 

We now summarize the results of the threshold analysis. To begin with, if the 
photon dispersion is unmodified and the electron parameter rj (for one helicity) 
is positive, then the electron group velocity Vg = 1 — {m? /2p^) + {rjp/M) + ■ • • 
exceeds the speed of light when 

Pth = {m^M/2r]f/^ ~ llTeV r]'^'^ . (C.21) 



This turns out to be the threshold energy for the vacuum Cerenkov process 
with emission of a zero energy photon, which we call the soft Cerenkov thresh- 
old. There is also the possibility of a hard Cerenkov threshold [53, 72]. For ex- 
ample, if the electron dispersion is unmodified and the photon parameter ^ is 
negative then at sufficiently high electron energy the emission of an energetic 
positive helicity photon is possible. This hard Cerenkov threshold occurs at 
Pth — (— 4m^M/^)^/^, and the emitted photon carries away half the incoming 
electron momentum. It turns out that the threshold is soft when both > 
and ^ > — 3?7, while it is hard when both ^ < —3ri and ^ < rj. The hard 
threshold in the general case is given by pth = (— 4m^M(^ + v)/i^ ~ vYY^^^ 
and the photon carries away a fraction — r))/2{$, + rj) of the incoming mo- 
mentum. In the general case at threshold, neither the incoming nor outgoing 
electron group velocity is equal to the photon phase or group velocity, so the 
hard Cerenkov effect cannot simply be interpreted as being due to faster than 
light motion of a charged particle. 



^■^ An extensive investigation of vacuum Cerenkov radiation was carried OTit in 
Ref. [110], focusing on the case of a Maxwell-Chern-Simons, dimension three LV 
operator for the photon. Unfortunately, that form of LV is different from what we 
study here, hence we cannot make direct use of those results. 
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C. 6 Fermion pair emission threshold 



The threshold for the process e~ e~e~e'^ involves only one LV parameter 
when the electrons all have the same helicity and the positron has the opposite 
helicity. In the threshold configuration this is allowed by angular momentum 
conservation, and the form of the spinors (14,17) shows that the amplitude is 
large as well, so rate will be high above threshold. It thus suffices to determine 
the threshold. 

Suppose the electron has positive helicity, and suppose that 7/+ > 0, so the 
reaction will proceed in this configuration. Let the incoming electron have 
momentum p, and let the positron have momentum zp. Since the electron 
dispersion relation then has positive curvature everywhere, as discussed in 
Appendix B the two final state electrons will have the same momentum, (1 — 
z)p/2. Writing out the conservation of energy in the threshold configuration, 
and keeping as usual just the first order terms in m? and ?7_|_, one finds the 
relation 

= — n — V *^-22 

r]+z[l - z) 



The right hand side is minimized when z = 1/2, yielding the threshold formula 
Pth = {%m^M/r^f'\ (C.23) 

which is a factor 16^/^ ~ 2.5 above the soft Cerenkov theshold (C.21). 

C.7 Helicity decay rate 



A variation on the Cerenkov effect that has received almost no attention in 
the literature is "helicity decay". If ?7_|_ and rj^ are unequal, say r/_ > 77+, then 
a negative helicity electron can decay into a positive helicity electron and a 
photon, even when the LV parameters do not permit the vacuum Cerenkov 
effect. In this process, the large R or small {0{m/E)) L component of a 
positive helicity electron is coupled to the small R or large L component of a 
negative helicity electron respectively. Such helicity decay has no threshold 
energy, so whether this process can be used to set a constraint is solely a 
matter of the decay rate. As we shall see however, the rate is maximum at an 
"effective threshold" energy below which it is very strongly suppressed, and 



Note that the relative size of L and R components of the fermion spinors is still 
controlled by the mass, as can be seen from Eq. 9. 
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above which it decreases as 1/E. Hence in practice there is a threshold for 
hehcity decay as welL 

The hehcity decay rate involves all three LV parameters ri± and ^. However, 
since the birefringence constraints limit ^ to a level far beyond where helicity 
decay is sensitive, we set ^ = for the purposes of evaluating the helicity 
decay rate here. 

We speciahze to the case of a negative helicity electron decaying into a positive 
helicity electron, the opposite case works identically. Again, the first step is 
to calculate the matrix element for the reaction, given by 

iM = ieM+(p')ea*7"M-(p) (C.24) 



where p^p' are the incoming and outgoing electron momenta, respectively. The 
calculation proceeds along the lines of photon decay, and the matrix element 
can be evaluated approximately as 




(C.25) 



No angular dependence occurs in this lowest order approximation for M since 
the helicity flip allows angular momentum to be conserved when all momenta 
are parallel. 

There are two different regimes for the size of the matrix element and rate, 
separated by the momentum pth ~ {m^M/ (rj^ — 77+))^/^. To see this, we must 
turn to the conservation equation. We will be interested in the amount of 
longitudinal phase space available to the reaction, i.e. the bounds p'^^ and p'^2 
in Eq. C.5. The energy conservation equation without any transverse momenta 
is 



where p' and k are the longitudinal components of the final electron and 
photon momenta. We introduce the variable z hj k = pz,p' = p{l — z). The 
photon energy is \k\, so energy conservation implies \z\ < 1. Since this is not a 
threshold situation, the final momenta can be anti-parallel, so negative values 
of z are permitted. In terms of z, Eq. C.26 becomes 

(C.27) 



m 



p_ 
.M 



{r^^-r^^{l-zf)-2{\z\-z) 
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li z < then the lowest order solution (using m^/p^ <g 1 and p/M -C 1) is 
z = —{rj- — rj+)p{4M)~^ . This is negative only if ?7_ > r/+ and is a solution 
for any p. Hence, there is no absolute threshold. In practice, p/M is so small 
that we can neglect the difference between this value of z and zero. We shall 
therefore set the lower bound on our rate integral for helicity decay to zero. 

Now consider ^ > 0. Eq. C.27 can in this case be rewritten as 



p3 



(C.28) 



At low energies {m^M/p^ ^ (77- — 77+)), the value of z that solves Eq. C.28 is 
also very small, z ^ (?7_ — r]+)p^ / {m^ M) = z. The corresponding bounds on 
the integral (C.5) arc therefore zero and z, so there is only a very small amount 
of phase space available for the reaction. Furthermore, since z is small, E' ^ E 
and (C.25) is doubly suppressed, once by the mass and once by AE — E — E'. 
Evaluating M explicitly at z — z gives 

3 

iM^ie{r,^-r,+)-^. (C.29) 



We can overestimate the rate at low energies by substituting (C.29) and the 
bounds on z into (C.5). It is an overestimate since at z — z the photon has its 
maximum possible energy, so the difference between E and E' is maximized, 
hence we are using the largest possible value of M in the allowed region of 
phase space. Upon integration the decay rate is at most 



If p » {m^M/ (rj- — r}+)y^^ then the situation is much different. The solution 
to the conservation equation (C.28) has z ~ 1 so the entire longitudinal phase 
space has opened up. Therefore E' ^ E and the only suppression is by the 
mass. The resulting decay rate in this case is approximately 



(C.31) 



which roughly coincides with Fi at the transition momentum — {rri^M/ (r)^ — 

V+)Y^'^ (as it must). Note that the decay rate actually decreases with increas- 
ing momentum. This is because the electrons become more and more chiral 
as momentum increases, thereby reducing their ability to flip helicity via the 
QED vertex. 
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For ?7_ — ?7+ of 0(1), ptr is approximately 10 TeV. Hence a 1 TeV negative 
helicity electron has a decay rate given by Fi which yields a lifetime of about 
one second. In contrast a 50 TeV electron with a decay rate of F2 has a 
lifetime of approximately 10~^ seconds, making a simple "threshold" analysis 
applicable for deriving constraints if 50 TeV negative helicity electrons are 
required to exist for the longer time scales relevant in a system such as the 
Crab nebula. 
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